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Abstract 



This article focuses on the two- flavor color superconducting phase at mod- 
erate baryon density. In order to simultaneously investigate the chiral phase 
transition and the color superconducting phase transition, the Nambu-Gorkov 
formalism is extended to treat the quark-antiquark and diquark condensates 
on an equal footing. The competition between the chiral condensate and the 
diquark condensate is analyzed. The cold dense charge neutral two-flavor 
quark system is investigated in detail. Under the local charge neutrality 
condition, the ground state of two-flavor quark matter is sensitive to the cou- 
pling strength in the diquark channel. When the diquark coupling strength is 
around the value obtained from the Fierz transformation or from fitting the 
vacuum bayron mass, the ground state of charge neutral two- flavor quark mat- 
ter is in a thermal stable gapless 2SC (g2SC) phase. The unusual properties 
at zero as well as nonzero temperatures and the chromomagnetic properties 
of the g2SC phase are reviewed. Under the global charge neutrality condition, 
assuming the surface tension is negligible, the mixed phase composed of the 
regular 2SC phase and normal quark matter is more favorable than the g2SC 
phase. A hybrid nonstrange neutron star is constructed. 
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I. INTRODUCTION 



A. QCD phase structure 

Quantum Chromodynamics (QCD) is regarded as the fundamental theory of quarks and 
gluons. Its ultimate goal is to explain all strong interaction experiments at all energies, high 
and low. QCD is an asymptotically free theory [1]. At very high energies, interaction forces 
become weak, thus perturbation calculations can be used. The perturbative QCD predictions 
have been extensively confirmed by experiments, while QCD in the non-perturbative regime 
is still a challenge to theorists. The fundamental quarks and gluons of QCD have not been 
seen as free particles, but are always confined within hadrons. It is still difficult to construct 
the hadrons in terms of nearly massless quarks and gluons. The observed baryon spectrum 
indicates that the (approximate) chiral symmetry is spontaneously broken in the vacuum. 
As a result, the eight pseudoscalar mesons tt, K and 77 are light pseudo-Nambu-Goldstone 
bosons, and the constituent quark obtains dynamical mass, which contributes to the baryon 
mass. Besides conventional mesons and baryons, QCD itself does not exclude the existence of 
the non- conventional states such as glueballs, hybrid mesons and multi-quark states [2] . Last 
year, pentaquarks, e.g., Q^iv^cPs) [3] with the mass M = 1540±10MeV and S {s^cPu) [4] 
with M — 1862±18MeV, were discovered in experiments. This has stirred a great interest of 
theorists to understand the structure of the pentaquarks [5-8] as well as other fundamental 
problems of QCD in the non-perturbative regime. 

Since 1970s, people have been interested in QCD at extreme conditions. It is expected 
that the chiral symmetry can be restored, and quarks and gluons will become deconfined at 
high temperatures and/or densities [9-12]. Results from lattice show that the quark- gluon 
plasma (QGP) does exist. For the system with zero net baryon density, the deconfinement 
and chiral symmetry restoration phase transitions happen at the same critical temperature 
[13]. At asymptotically high temperatures, e.g., during the first microseconds of the "Big 
Bang" , the many-body system composed of quarks and gluons can be regarded as an ideal 
Fermi and Boson gas. It is believed that the "httle Bang" can be produced in RHIC and 
LHC. Recently, it was shown that the new state of matter produced in RHIC is far away 
from the asymptotically hot QGP, but in a strongly coupled regime. This state is called 
strongly coupled quark-gluon plasma (sQGP). In this strong coupling system, the meson 
bound states still play an important role [14,15]. For most recent reviews about QGP, e.g., 
see Ref. [16-19]. 

Studying QCD at finite baryon density is the traditional subject of nuclear physics. The 
behaviour of QCD at finite baryon density and low temperature is central for astrophysics to 
understand the structure of compact stars, and conditions near the core of collapsing stars 
(supernovae, hypernovae). Cold nuclear matter, such as in the interior of a Pb nucleus, 
is at T = and hb — = 940MeV. Emerging from this point, there is a first-order 
nuclear liquid-gas phase transition, which terminates in a critical endpoint at a temperature 
~ lOMeV [20]. If one squeezes matter further and further, nuleons will overlap. Quarks and 
gluons in one nucleon can feel quarks and gluons in other nucleons. Eventually, deconfine- 
ment phase transition will happen. Unfortunately, at the moment, lattice QCD is facing 
the "sign problem" at nonzero net baryon densities. Our understanding at finite baryon 
densities has to rely on effective QCD models. Phenomenological models indicated that, at 
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nonzero baryon density, the QGP phase and the hadron gas are separated by a critical hne 
of roughly a constant energy density — IGeV/fm'^ [21]. 

In the case of asymptotically high baryon density, the system is a color superconducor. 
This was proposed 25 years ago by Prautschi [22] and Barrois [23]. Reaching this conclusion 
does not need any other knowledge, if one knows that one-gluon exchange between two 
quarks is attractive in the color antitriplet channel, and if one is also familiar with the theory 
of superconductivity, i.e., the Bardeen, Cooper, and Schrieffer (BCS) theory [24]. Based on 
the BCS scenario, if there is a weak attractive interaction in a cold Fermi sea, the system 
is unstable with respect to the formation of particle-particle Cooper-pair condensate in the 
momentum space. Detailed numerical calculations of color superconducting gaps were firstly 
carried out by Bailin and Love [25]. They concluded that the one-gluon exchange induces 
gaps on the order of 1 MeV at several times of nuclear matter density. This small gap 
has little effect on cold dense quark matter, thus the investigation of cold quark matter lay 
dormant for several decades. It was only revived recently when it was found that the color 
superconducting gap can be of the order of 100 MeV [26], which is two orders larger than 
early perturbative estimates in Ref. [25]. For this reason, the topic of color superconductivity 
stirred a lot of interest in recent years. For review articles on the color superconductivity, 
see for example, Refs. [16,27-33]. 

In this paper, I will focus on the color superconducting phase structure at intermediate 
baryon density regime based on my own experience. It is worth to mention that there is 
another good review article Ref. [32] in this regime from another point of view. The outline 
of this article is as follows: A brief overview of color superconducting phases is given in Sec. 
I B. In Sec. II, the generalized Nambu-Gorkov formalism is introduced to treat the chiral and 
diquark condensates on an equal footing, and the competition between the chiral and color 
superconducting phase transitions is investigated. Sec. Ill focuses on homogeneous neutral 
two-flavor quark matter, the gapless 2SC (g2SC) phase and its properties are reviewed. In 
Sec. IV, a neutral mixed phase composed of the 2SC phase and normal quark matter is 
discussed. At last, a brief outlook is given in Sec. V. 



B. Color superconducting phases 

Let us start with the system of free fermion gas. Fermions obey the Pauli exclusion 
principle, which means no two identical fermions can occupy the same quantum state. The 
energy distribution for fermions (with mass m) has the form of 

here Ep — -^p^ _|_ ^ jg ^^le chemical potential and T is the temperature. At zero 
temperature, f{Ep) — 9{iJ, — Ep). The ground state of the free fermion gas is a filled 
Fermi sea, i.e., all states with the momenta less than the Fermi momentum pp — \J p?' — 

are occupied, and the states with the momenta greater than the Fermi momentum pp are 
empty. Adding or removing a single fermion costs no free energy at the Fermi surface. 

For the degenerate Fermi gas, the only relevant fermion degrees of freedom are those 
near the Fermi surface. Considering two fermions near the Fermi surface, if there is a net 
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attraction between them, it turns out that they can form a bound state, i.e.. Cooper pair 
[34]. The binding energy of the Cooper pair A{K) {K the total momentum of the pair), is 
very sensitive to K, being a maxium where K = 0. There is an infinite degeneracy among 
pairs of fermions with equal and opposite momenta at the Fermi surface. Because Cooper 
pairs are composite bosons, they will occupy the same lowest energy quantum state at zero 
temperature and produce a Bosc-Einstein condensation. Thus the ground state of the Fermi 
system with a weak attractive interaction is a complicated coherent state of particle-particle 
Cooper pairs near the Fermi surface [24] . Exciting a quasiparticle and a hole which interact 
with the condensate requires at least the energy of 2 A. 

In QED case in condensed matter, the interaction between two electrons by exchanging 
a photon is repulsive. The attractive interaction to form electron-electron Cooper pairs is by 
exchanging a phonon, which is a collective excitation of the positive ion background. The 
Cooper pairing of the electrons breaks the electromagnetic gauge symmetry, and the photon 
obtains an effective mass. This indicates the Meissner effect [35], i.e., a superconductor 
expels the magnetic fields. 

In QCD case at asymptotically high baryon density, the dominant interaction between 
two quarks is due to the one-gluon exchange. This naturally provides an attractive interac- 
tion between two quarks. The scattering amplitude for single-gluon exchange in an SU {Nc) 
gauge theory is proportional to 

{Ta)ki{Ta)lj = - SikSjl) + {SjkSii + SikSji). (1.2) 

Where is the generator of the gauge group, and i,j and k,l are the fundamental colors 
of the two quarks in the incoming and outgoing channels, respectively. Under the exchange 
of the color indices of either the incoming or the outgoing quarks, the first term is antisym- 
metric, while the second term is symmetric. For Nc — 3, Eq. (1.2) represents that the tensor 
product of two fundamental colors decomposes into an (antisymmetric) color antitriplet and 
a (symmetric) color sextet, 

[3]^ ® [3]^ = [3]^ © [6]^. (1.3) 

In Eq. (1.2), the minus sign in front of the antisymmetric contribution indicates that the 
interaction in this antitriplet channel is attractive, while the interaction in the symmetric 
sextet channel is repulsive. 

For cold dense quark matter, the attractive interaction in the color antitriplet channel 
induces the condensate of the quark-quark Cooper pairs, and the ground state is called the 
"color superconductivity" . Since the diquark cannot be color singlet, the diquark condensate 
breaks the local color SU{3)c symmetry, and the gauge bosons connected with the broken 
generators obtain masses. Comparing with the Higgs mechanism of dynamical gauge sym- 
metry breaking in the Standard Model, here the diquark Cooper pair can be regarded as a 
composite Higgs particle. The calculation of the energy gap and the critical temperature 
from the first principles has been derived systematically in Refs. [36-43]. 

In reality, we are more interested in cold dense quark matter at moderate baryon density 
regime, i.e., Hq ~ 500MeV, which may exist in the interior of neutron stars. It is likely that 
cold dense quark droplet might be created in the laboratory through heavy ion collisions 



4 



in GSI-SPS energy scale. At these densities, an extrapolation of the asymptotic arguments 
becomes unreliable, we have to rely on effective models. Calculations in the framework of 
pointlike four-fermion interactions based on the instanton vertex [26,44-46], as well as in the 
Nambu-Jona-Lasinio (NJL) model [47-51] show that color superconductivity does occur at 
moderate densities, and the magnitude of diquark gap is around 100 MeV. 

Even though the antisymmetry in the attractive channel signifies that only quarks with 
different colors can form Cooper pairs, color superconductivity has very rich phase structure 
because of its flavor, spin and other degrees of freedom. In the following, I list some of the 
known color superconducting phases. 

The 2SC phase 

Firstly we consider a system with only masslcss u and d quarks, assuming that the 
strange quark is much heavier than the up and down quarks. The color superconducting 
phase with only two flavors is normally called the 2SC phase. 

Renormalization group arguments [36,52,53] suggest that possible quark pairs always con- 
dense in the s— wave. This means that the spin wave function of the pair is anti-symmetric. 
Since the diquark condenses in the color antitriplet 3c channel, the color wave function of 
the pair is also anti-symmetric. The Pauli principle requires that the total wave function 
of the Cooper pair has to be antisymmetric under the exchange of the two quarks forming 
the pair. Thus the flavor wave function has to be anti-symmetric, too. This determines the 
structure of the order parameter 

= A6,,6"^^ (1.4) 

where color indices a,f3 E {r,g,b) and flavor indices i,j G {u,d). From the order parameter 
Eq. (1.4), we can see that the condensate picks a color direction (here the blue direction, 
which is arbitrarily selected). The ground state is invariant under an SU{2)c subgroup of 
the color rotations that mixes the red and green colors, but the blue quarks are singled out 
as different. Thus the color SU{3)c is broken down to its subgroup SU{2)c, and five of the 
gluons obtain masses, which indicates the Mcissner effect [54]. 

In the 2SC phase, the Cooper pairs are ud — du singlets and the global ffavor symmetry 
SU {2)l®SU {2)r is intact, i.e., the chiral symmetry is not broken. There is also an unbroken 
global symmetry which plays the role of C/(1)b. Thus no global symmetry are broken in the 
2SC phase. 

The CFL phase 

In the case when the chemical potential is much larger than the strange quark mass, 
we can assume m„ = = = 0, and there are three degenerate massless flavors in the 
system. The spin-0 order parameter should be color and flavor anti-symmetric, which has 
the form of 

Ay = A5^6,,,6"^^ (1.5) 
I 

where color indices a,(3 & {r,g,b) and flavor indices i,j G {u,d,s). Writing ^/Cij/e"'^^ = 
SfSj — Sj'S^, we can see that the order parameter 
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Ay = A(5r5f - 5j^5f ) (1.6) 

describes the color-flavor locked (CFL) phase proposed in Ref. [55]. Many other different 
treatments [56-58] agreed that a condensate of the form (1.5) is the dominant condensate 
in three-flavor QCD. 

In the CFL phase, all quark colors and flavors participate in the pairing. The color gauge 
group is completely broken, and all eight gluons become massive [55,59], which ensures that 
there are no infrared divergences associated with gluon propagators. Electromagnetism is 
no longer a separate symmetry, but corresponds to gauging one of the flavor generators. A 
rotated electromagnetism ("Q") remains unbroken. 

Two global symmetries, the chiral symmetry and the baryon number, are broken in the 
CFL phase, too. In zero-density QCD, the spontaneous breaking of chiral symmetry is due 
to the condensation of left-handed quarks with right-handed quarks. Here, at high baryon 
density, the chiral symmetry breaking occurs due to a rather different mechanism: locking 
of the flavor rotations to color. In the CFL phase, there is only pairing of left-handed quarks 
with left-handed quarks, and right-handed quarks with right-handed quarks, i.e., 

< rLi^% >= - < rRi^% > ■ (1-7) 

Where L, i? indicate left- and right-handed, respectively, a,/3 are color indices and i,j are 
flavor indices. A gauge invariant form [27,60,61] 

< ^£.<^LV^W >~< X i'lai'Rp >~ A^e^.-^e^'- (1.8) 

captures the chiral symmetry breaking. The spectrum of excitations in the CFL phase 
contains an octet of Goldstone bosons associated with the chiral symmetry breaking. This 
looks remarkably hke those at low density. In the excitation spectrum of the CFL phase, 
there is another singlet U{1) Goldstone boson related to the baryon number symmetry 
breaking, which can be described using the order parameter 

< udsuds >~< AA > . (1.9) 

In QCD with three degenerate light flavors, the spectrum in the CFL phase looks similar to 
that in the hyper-nuclear phase at low-density. It is suggested that the low density hyper- 
nuclear phase and the high density quark phase might be continuously connected [62]. 

Spin-l color superconductivity 

In the case of only one-flavor quark system, due to the antisymmetry in the color space, 
the Pauli principle requires that the Cooper pair has to occur in a symmetric spin channel. 
Therefore, in the simplest case, the Cooper pairs carry total spin one. Spin-l color super- 
conductivity was firstly studied in Ref. [25], for more recent and detailed discussions about 
the spin-l gap, its critical temperature and Meissner effect, see Refs. [38,39,63-66]. For a 
review, see Ref. [67]. 

Pairing with mismatch: LOFF, CFL-K, g2SC and gCFL 

To form the Cooper pair, the ideal case is when the two pairing quarks have the same 
Fermi momenta, i.e., pF,i — Pfj with p^j = \ l^'pi ~ "^f' ^^^^ ™ ideal 2SC, CFL, and 
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spin-1 color superconducting phases. However, in reality, the nonzero strange quark mass 
or the requirement of charge neutrality induces a mismatch between the Fermi momenta of 
the two pairing quarks. When the mismatch is very small, it has little effect on the Cooper 
pairing. While if the mismatch is very large, the Cooper pair will be destroyed. The most 
interesting situation happens when the mismatch is neither very small nor very large. 

LOFF: In the regime just on the edge of decoupling of the two pairing quarks (due 
to the nonzero strange quark mass for the qs Cooper pair with q G {u, d) or the chemical 
potential difference for the ud Cooper pair), a "LOFF" (Larkin-Ovchinnikov-Fulde-Ferrell) 
state may be formed. The LOFF state was firstly investigated in the context of electron 
superconductivity in the presence of magnetic impurities [68,69]. It was found that near the 
unpairing transition, it is favorable to form a state in which the Cooper pairs have nonzero 
momentum. This is favored because it gives rise to a regime of phase space where each of 
the two quarks in a pair can be close to its Fermi surface, and such pairs can be created 
at low cost in free energy. This sort of condensates spontaneously break translational and 
rotational invariance, leading to gaps which vary periodically in a crystalline pattern. The 
crystalline color superconductivity has been investigated in a series of papers, e.g., see Refs. 
[70-77]. 

CFL-K: The strange quark mass induces an effective chemical potential /ig — 
ml/{2pp), and the effects of the strange quark mass can be quite dramatic. In the CFL 
phase, the K'^ and K'^ modes may be unstable for large values of the strange quark mass 
to form a kaon condensation [78-81]. In the framework of effective theory [60,61,82-84], the 
masses of the Goldstone bosons can be determined as 



m: 



T- 



2PF 



+ 
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with A = 3A^/(47r^) [61,85]. It was found that the kaon masses are substantially affected 
by the strange quark mass, the masses of and are pushed up while and are 
lowered. As a result, the and K'^ become massless if ms\crit = 3.03 m^/^A'^^^. For larger 
values of the kaon modes are unstable, signaling the formation of a kaon condensate. 
For review of the kaon condensate in the CFL phase, see Ref. [31]. Recently, it was found 
that in the CFL phase, there also may exist r] condensate [86]. 

g2SC and gCFL: When the /5-equilibrium and the charge neutrality condition are 
required for the two-flavor quark system, the Fermi surfaces of the pairing u quark and d 
quark differ by He, here /ig is the chemical potential for electrons. It was found that when 
the gap parameter A < /Xe/2, the system will be in a new ground state called the gapless 
2SC (g2SC) phase [87]. The g2SC phase has very unusual temperature properties [88] and 
chromomagnetic properties [89]. This phase will be introduced in more detail in Sec. 111. 

Similarly, for a charge neutral 3-fiavor system with a nonzero strange quark mass m^, 
with increasing m^, the CFL phase transfers to a new gapless CFL (gCFL) phase when 
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m^/lJ, ~ 2 A [90]. The finite temperature property of tlie cliarge neutral tliree-fiavor quark 
matter was investigated in Ref. [91-93]. Recently, it was shown that the kaon condensate 
shifts the critical strange quark mass to higher values for the appearance of the gCFL phase 
[94]. 

II. COMPETITION BETWEEN THE CHIRAL AND DIQUARK CONDENSATES 

In this section, we are going to investigate the competition between the chiral and diquark 
condensates, and the content is based on Refs. [51,97]. 

In the idealized case at asymptotically high baryon densities, the color superconductivity 
with two massless flavors and the color-flavor-locking (CFL) phase with three degenerate 
massless quarks have been widely discussed from first principle QCD calculations. Usually, 
the diagrammatic methods are used at the asymptotic densities. The Green-function of the 
eight-component field and the gap equation were discussed in details in [38,95]. Neither 
current quark mass nor chiral condensate are necessary to be considered because they can 
be neglected compared to the very high Fermi surface. At less-than-asymptotic densities, 
the corrections of nonzero quark mass to the pure CFL phase can be treated perturbatively 
by expanding the current quark mass around the chiral limit [85,96]. 

For physical applications, we are more interested in the moderate baryon density regime, 
which may be related to the neutron stars and, in very optimistic cases even to heavy-ion 
collisions. To work out the phase structure from the hadron phase to the color super- 
conducting phase, one should deal with the chiral condensate and the diquark condensate 
simultaneously. Because the chiral condensation contributes a dynamical quark mass, it is 
not reasonable any more to treat the quark mass term perturbatively in this density regime. 

Usually, at moderate baryon density regime, effective models such as the instanton model, 
as well as the Nambu-Jona-Lasinio (NJL) model, are used. The model parameters are fixed 
in the QCD vacuum. In this regime, the usual way is to use the variational methods, which 
are enough for constructing the thermodynamic potential and the gap equation to investigate 
the ground state. In order to investigate the diquark excitations, and to calculate the cross- 
section of physical process in the color superconducting phase, it is very helpful to use the 
Green-function approach. 

One of our main aims in this section is to apply the Green-function approach in the 
moderate baryon density regime. By using the energy projectors for massive quark, we can 
deal with the chiral and diquark condensates on an equal footing, and we will deduce the 
Nambu-Gorkov massive quark propagator [97] as well as the thermodynamic potential [51]. 

In the normal phase, the quarks of different colors are degenerate. While in the color 
broken phase, it is natural to assume that the quarks involved in the diquark condensate 
are different from the free quarks. It is difficult to obtain the mass expression for the quarks 
participating in the diquark condensate, because the particles and holes mix with each other, 
and the elementary excitations are quasi-particles and quasi-holes near the Fermi surface. 
The difference between quarks in different colors has been reflected by their propagators, 
and the difference can be read through calculating the quarks chiral condensate. In the 
chiral limit, the chiral condensate disappears entirely in the color superconducting phase, 
and it is not possible to investigate the influence of the color breaking on quarks in different 
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colors, so we keep the small current quark mass in this section. 

In the moderate baryon density regime, people are interested in whether there exists 
a regime where both the chiral symmetry and color symmetry are broken [45-47,98-102]. 
Note that the presence of a small current quark mass induces that the chiral symmetry 
only restores partially and there always exists a small chiral condensate in the color super- 
conducting phase. This phenomenon had been called coexistence of the chiral and diquark 
condensates in Ref. [45]. In the coexistence regime resulted by current quark mass, the chiral 
condensate is small comparing with the diquark condensate, and the role of the chiral con- 
densate can be prescribed by the Anderson theorem [102], i.e., in this phase, the contribution 
of the chiral condensate to thermodynamic quantities becomes strongly suppressed, and one 
can calculate the diquark condensate neglecting the influence of the chiral condensate. 

In order to differ from the coexistence of the diquark condensate and the small chiral 
condensate in the chiral symmetry restoration phase, the coexistence of the diquark conden- 
sate and the dynamical chiral condensate in the chiral symmetry breaking phase is called 
the double broken phase. The existence of the double broken phase depends on the cou- 
pling constants Gs and Go in the quark- ant iquark and diquark channels [47,101]. In the 
case of small ratio of Gd/Gs < 1, the calculations in the instanton model [45], NJL model 
[47] and random matrix model [101] show that there is no double broken phase. While for 
Gd/Gs > 1 the calculations in the the random matrix model [101] and the NJL model 
[47] show that there exists a narrow regime where both chiral (dynamical) symmetry and 
color symmetry are broken, and the chiral and diquark condensate coexist. The larger value 
of Gd/Gs is, the wider regime of the double broken phase has been found in the random 
matrix model [101]. 

In this section, based on Ref. [51], we explain the existence of the double broken phase 
by analyzing the influence of the diquark condensate on the Fermi surface of the constituent 
quark. In the mean-field approximation of the NJL model, the thermal system of the 
constituent quarks is nearly an ideal Fermi gas, and there is a sharp Fermi surface for 
the constituent quark. When a diquark condensate is formed, the Cooper pair smoothes the 
Fermi surface, and this induces the chiral symmetry restoring at a smaller chemical potential. 
The stronger the coupling constant in the diquark channel, the larger the diquark gap will 
be, and the smaller critical chemical potential will be for the color superconductivity phase 
transition. However, it is found in Ref. [51] that, if Gd/Gs > 1, the diquark condensate 
starts to appear even in the vacuum. This may suggest that in the NJL model, the diquark 
coupling strength Gd should not be larger than the quark-antiquark coupling strength Gs- 
It is also possible that the appearance of the diquark condensate in vacuum is an artificial 
result due to the lack of deconfinement in the NJL model. 

A. The extended NJL model 

To describe quark matter in the intermediate baryon density regime, we use the extended 
NJL model. The choice of the NJL model [103] is motivated by the fact that this model 
displays the same symmetries as QCD and that it describes well the spontaneous breakdown 
of chiral symmetry in the vacuum and its restoration at high temperature and density. For 
the review of the NJL model, please see [104]. The model we used in this paper is an 
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extended version of the two-flavor NJL model inchiding interactions in the color singlet 
quark- ant iquark channel as well as in the color anti-triplet diquark channel, which is not 
directly extended from the NJL model, but from the QCD Lagrangian [105-107]. 

The importance of color 3 diquark degree of freedom is related to the fact that one can 
construct a color-singlet nucleon current based on it. Because the gluon exchange between 
two quarks in the color 3 channel is attractive, one can view a color singlet baryon as a 
quark-diquark bound state. And experimental data from ep collisions indicate the existence 
of this quark-diquark component in nucleons [108]. The discovery of the pentaquark also 
reminds us the importance of the diquark degree of freedom [6,7,109]. 

The first attempt to investigate the diquark properties in the NJL model was taken in 
Ref. [110]. Starting from an NJL model for scalar, pseudoscalar, vector and axial- vector 
interactions of the {qq) x (qq) type and Fierz-transforming away the vector and axial- vector 
interactions, the scalar and pseudoscalar mesons, and diquarks can be obtained. However, 
this method could not get a consistent treatment of vector and axial-vector particles. 

The extended NJL model we use is derived from the QCD Lagrangian [105,106]. Inte- 
grating out the gluon degrees of freedom from the QCD Lagrangian, and performing a local 
approximation for the (nonperturbative) gluon propagator, one obtains a contact current- 
current interaction. By using a special Fierz-rearrangement [111] (see Appendix A), one can 
completely decompose the two-quark-current interaction term into "attractive" color singlet 
(qq) and color antitriplet (qq) channels. In this way, a complete simultaneous description 
of scalar, pseudo-scalar, vector, and axial-vector mesons and diquarks is possible, thus the 
extended NJL model including {qq) x {qq) interactions is completed by a corresponding 
{qq) X {qq) interaction part. 

It is worth to mention that, in the mean- field approximation, there is a famous am- 
biguity connected with performing the Fierz transformation for the pointlike four-fermion 
interaction. For more detail discussion about the Fierz ambiguity, see Ref. [113]. Here, in 
this paper, we only consider the scalar, pseudoscalar mesons and scalar diquark, and the 
Lagrangian density has the form 

C = q{ird^ - mo)q + Gs[{qq)' + {qi^^rq)'] + GD[{ifee'^,q){iqse'^,q'')]. (2.1) 

Where q^ — Cq^, q^ — q^C are charge- conjugate spinors, C — ^7^7° is the charge con- 
jugation matrix (the superscript T denotes the transposition operation). For more details 
about the charge conjugate, see Appendix B. mo is the current quark mass, the quark field 
q = q^a with i = 1,2 and a = 1,2,3 is a fiavor doublet and color triplet, as well as a four- 
component Dirac spinor, r = (t^,t^,t^) are Pauli matrices in the flavor space, where is 
antisymmetric, and {eY'' = e^'', {e^)°'^ = e'*^^ are totally antisymmetric tensors in the flavor 
and color spaces. 

Though, from the Fierz transformation, the quark-antiquark coupling constant Gs and 
the diquark coupling constant Gd have the relation Gd = Gs Nc/{2Nc — 2), in this article, 
Gs and Go are treated independently. The former is responsible for the meson excitations, 
and the latter for the diquark excitations, which in principle can be determined by fltting 
meson and baryon properties in the vacuum. In Ref. [106], the ratio of coupling constants 
Gd/Gs — 2.26/3 was obtained by fitting the scalar diquark mass of 600 MeV to get a 
realistic vacuum baryon mass. It is noticed that this ratio is quite close to the value 3/4 
obtained from Fierz transformation. 
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The attractive interaction in different cliannels in tlic Lagrangian gives rise to a very rich 
structure of the phase diagram. At zero temperature and density, the attractive interaction 
in the color singlet channel is responsible for the appearance of a quark-antiquark condensate 
and for the spontaneous breakdown of the chiral symmetry, and the interaction in the qq 
channel binds quarks into diquarks (and baryons), but is not strong enough to induce diquark 
condensation. As the density increases, Paiili blocking suppresses the qq interaction, while 
the attractive interaction in the color anti-triplet diquark channel will induce the quark- 
quark condensate around the Fermi surface which can be identified as a superconducting 
phase. 

After bosonization [105,106,112], one obtains the linearized version of the model, i.e., in 
the mean-field approximation. 



C = q{iYd^ - mo)q - q{a + i'y^T7r)q - ^A*\ifee%q) - ^A^iqee^^^q^) 

^2 ^ ^2 ^*b^b 



(2.2) 



with the bosonic fields 

iq^ee^^^q, A*^ iqee^^^q^ , a qq, n iq^^rq. (2.3) 



Clearly, the a and tt fields are color singlets, and the diquark fields A** and A**" are color 
antitriplet and (isoscalar) singlet under the chiral SU (2) l x SU (2)r group, a ^ Q and A^ 
indicate that the chiral symmetry and the color symmetry are spontaneously broken. Here 
it has been regarded that only the red and green quarks participating in the condensate, 
while the blue quarks do not. The real ground state at any T, will be determined by the 
minimum of the thermodynamic potential. 

The partition function of the grand canonical ensemble can be evaluated by using the 
standard method [114,115], 

Z = N' J [dq][dq]exp{ dr J d^x {£ + i^q^oq)}, (2.4) 

where fi is the chemical potential, and P — l/T is the inverse of the temperature T. 

In the mean field approximation, we can write the partition function as a product of 
three parts, 

Z — ZconstZbZr^g- (2-5) 

The constant part is 

Zconst = A^'exp{- dr j d^x + ^}. (2.6) 
The contribution from free blue quarks has the form 

= y [(igb][(igfe]exp{ J J d^x [^qb{iY9f, - m + ix^o)qb 

+\^{il'd^.-m-lil^)q^]}. (2.7) 
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For the quarks with red and green color Q — gr^.g participating in the quark condensate, 
their contribution is 



■7 

^r,g 



[dQ][dQ]exp{J^ dr J d^^ [^Qit^'^d^ - m + i^^o)Q + 

lO'^iiYd, -m- ix^o)Q'' + ^QA-g^ + ^g^A+g]}. (2.8) 
Here we have introduced the constituent quark mass 

m — mo + (7, (2.9) 

and defined A='= as 

A- = -iAeeS, A+ = -i/\*ee% (2.10) 
with the relation A+ = -f^{/Sr)'^-f^ . 

B. Nambu-Gorkov propagator with chiral and diquark condensates 

In the case of finite chemical potential, it is more convenient to use the Nambu-Gorkov 
formalism. Introducing the 8-component spinors for the blue quarks and the quarks colored 
with red and green, respectively 



*f=( ), *^>=(?;. (2.11) 



and using the Fourier transformation in the momentum space, 

?(^) = E E ^^-^^'^"^-^■^^^(p), (2.13) 

^ n p 

where V is the volume of the thermal system, we can re- write the partition function Eqs. (2.7) 
and (2.8) in the momentum space as 

n,p 

= Det^/^(/3Go '), (2.14) 

and 

n,p 

= Det^/'(/3G-^). (2.15) 
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Where the determinantal operation Det is to be carried out over the Dirac, color, flavor 
and the momentum- frequency space. In Eqs. (2.14) and (2.15), we have defined the quark 
propagator in the normal phase 

with 

[G±]-i = (po±//)7o-7-p-m (2.17) 
and the quark propagator in the color broken phase 

G-i^ ( [G^l ' ^" y (2.18) 
The Nambu-Gorkov propagator G(P) is determined from solving 1 = G~-'^G, 

G = ( S+ ^- ) , (2.19) 

with the components 

S± = -G^ A± G^ = -G^ A± G^. (2.20) 

Here all components depend on the 4-momentum P. 

In the case of the mass term m — 0, the Nambu-Gorkov quark propogator has a simple 
form which could be derived from the energy projectors for massless particles [38,95]. If 
there is a small mass term, the quark propagator can be expanded perturbatively around 
m = 0, but the form is very complicated [96]. In our case the quark mass term cannot 
be treated perturbatively, we have to find a general way to deal with the massive quark 
propagator. 

Fortunately, we can evaluate a simple form for the massive quark propagator by using 
the energy projectors for massive particles. The energy projectors onto states of positive 
and negative energy for free massive particles are defined as 

A.(p) = i(l±35(l^), (2.21) 

where the quark energy Ep — ^/p^Tm^. Under the transformation of 70 and 75, we can 
get another two energy projectors A±, 

A±(P) = 2^-^^ E ' ^ 
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which satisfy 

7oA±(p)7o = At(p), 75A±(p)75 = A±(p). (2.23) 
The normal quark propagator elements can be re-written as 



.± 7oA+ , 7oA_ 



with — Ep± fj,. The propagator has four poles, i.e., 

Po^±E;, Po^tE;, (2.25) 

where the former two correspond to the excitation energies of particles and holes, and the 
latter two are for antiparticles and antiholes, respectively. 

The full quark propagator G, i.e., (2.19) can be evaluated [97]. The normal quark 
propagator has the form 



= r27oA+ + ^ £27oA-) {5^p - 5^,8p,)5'^ , (2.26) 



and the abnormal quark propagator reads 

Po -Et pI- El 

withE^^ = Ef + ^^. This propagator is similar to the massless propagator derived in [95]. 
The four poles of the Nambu-Gorkov propagator, i.e., 

Po^^El, Po^tEX, (2.28) 

correspond to the excitation energies of quasi-particles (quasi-holes) and quasi-antiparticle 
(quasi- antiholes) in the color broken phase. These quasi-particles are superpositions of 
particles and holes, and are called "Bogoliubons" . In the normal phase, exciting a pair of a 
particle and a hole on the Fermi surface does not need energy, while in the superconducting 
phase, exciting a quasi-particle and a hole needs at least the energy 2A at Ep — /i. 



C. Gap equations 

With the Nambu-Gorkov quark propagator Eqs. (2.26) and (2.27), one can investigate 
the meson excitations [116], and physical processes related to the transport properties. The 
chiral and diquark condensates can also be generally expressed from the quasi-particle prop- 
agator. The diquark condensate has the form 

< f^^q >= -iT J2 I (^^^^[^"^s]. (2.29) 
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From general consideration, there should be eight Dirac components of the diquark conden- 
sates [25,95,117,118]. In the case of the NJL type model, the diquark condensates related 
to momentum vanish, and there is only one independent O"*" diquark gap with Dirac struc- 
ture r = 75 for massless quark, and there exists another 0"*" diquark condensate with Dirac 
structure F = 7075 at nonzero quark mass. In this paper, we assume the contribution of the 
diquark condensate with F = 7075 is small, and only consider the diquark condensate with 

r = 75- 

Performing the Matsubara frequency summation and taking the limit T — >^ 0, we get the 
diquark condensate at finite chemical potential 

<,-«>=-2A^.JV,/|E^l^.^l. (2.30) 

Using the relation between the diquark gap A and the diquark condensate < qP^^Q. >, i-e., 

A = -2Gd < gSsg >, (2.31) 
the gap equation for the diquark condensate in the limit of T ^ can be written as 

The quark mass has the expression 

m = niQ — 2Gs < qq >, 
<qq> = 2<qrqr> + < Mb > ■ (2.33) 

Where the chiral condensate for the free blue quarks is 

< Mb >= -iTj2 J ^M^oi (2-34) 
and the chiral condensate for the red and green quarks reads 

< qrqr >= -iTY^ I -^tr[G+]. (2.35) 

In the limit T — > 0, the gap equation for the quark mass has the form 

m = mo - SmNfGs / [^(/^ ' Ep) - 1 + 2« - n;)] . (2.36) 

Where 

El 
^1 



nt^^ilT^) (2.37) 



are the occupation numbers for quasi-particles and quasi-antiparticles at T = 0. Correspond- 
ingly, 1 — Up are the occupation numbers of quasi-holes and quasi- antiholes, respectively. 
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The difference between the quarks participating in the diquark condensate and the free 
blue quark can be read from their chiral condensates, 

5 =< Mb - < mr >^'\ (2.38) 

where 5 has the dimension of energy. In the case of chiral limit, the quark mass m decreases 
to zero in the color superconducting phase, and the influence of the diquark condensate on 
quarks in different colors vanishes. 



D. Derivation of the thermodynamic potential 

Now let us derive the thermodynamic potential. For the blue quarks which do not 
participate in the diquark condensate, from Eq. (2.14), we have 

InZ,, = lln{Det(/3[Go]~') = ^ln[Det(/5[G+]-i)Det(/5[Go ]"')]• (2-39) 
Using the Dirac matrix, we first perform the determinant in the Dirac space, 
Det (3 [G-^]-^ = Det (3 [(po + //)7o - 7 • P - 

= Det/3f (^° + ^)-"^ , 

\^ -cr-p -(po + /^)-m 

= [{p,^^f-Ell (2.40) 

and in a similar way, we get 

Det /? [Go = [(Po - - Ell (2.41) 

After performing the determinant in the Dirac space, we have 

Det /3 [G+]-^Det /? [Go]"^ = /3>^ - {E, + fxf] f3'\pl - {E, - fif]. (2.42) 

Considering the determinant in the flavor, color, spin spaces and momentum-frequency 
space, we get the standard expression 

InZ,, = NjY, EO^(/^'t^o - {E, + ^,f]) + ln(/3>g - {E, - /.)^])}, (2.43) 

n p 

remembering that the color space for the blue quark is one-dimensional. 

It is more complicated to evaluate the thermodynamic potential for the quarks partici- 
pating in the diquark condensate. Prom Eq. (2.15), we have 

InZ,^,^ = ^lnDet(/3G-i). (2.44) 
For a 2 X 2 matrix with elements A,B,C and D, we have the identity 

( G ^ ) " Det(-GB + CAC-^D) Det(-SG + BDE'^A). (2.45) 
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Replacing A,B,C and D with the corresponding elements of G ^, we have 

Det(/?G-^) = P^DeWi = /^^Detf-A+A" + A+[G+]-^[A+]-i[Go ]"!] 

= P^DetD^ = /52Det[-A+A- + [G^ ]-i[A-]-i[G+]-iA-]. (2.46) 

Using the energy projectors A±, we can work out Di and D2 as 

L>i = A' + 75[7o(po - E;)A+ + 7o(po + £^+)A-]75[7o(po - E+)A+ + 7o(po + E-)A_] 

= -[(Po - Ef - A^)A_ + (p^ - Ef - A^)A+], 

^2 = -[{pI - Ef - A2)A+ + {pI - Ef - A2)A_]. (2.47) 

Using the properties of the energy projectors, one can get 

D,D2 = \pl-Ef]\pl-Ef]. (2.48) 

With the above equations, Eq. (2.44) can be expressed as 

InZ,,,^ = ^ln[Det/3G-^] = -^TrXnlp^D.^^D^] 

= \{TM(3'{pI - Ef)] + TTln^pl - EX')]} 

- 2A^/ E J2i^^if^"(Po - ^a)] + H(3\pI - Ef)]}. (2.49) 



n p 



The frequency summation of the free-energy 

\nZf = J2H(^\pl-Ei)] (2.50) 

n 

can always be obtained by performing the frequency summation of the propagator I /{pi — 
Ep). Differentiate Eq. (2.50) with respect to Ep-. 

where f{x) = l/{e^^ + 1) is the usual Fermi-Dirac distribution function. Then integrating 
with respect to Ep, one can get the free-energy 

InZf = p[Ep + 2Tln(l + e'^^^)]. (2.52) 

With the help of the above expression, and replacing 



one can have 
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InZ,, = Nj(5V j ^^[E; + 2Tln(l + e'^^) + E; + 2Tln(l + e"'^^^")], 
InZ,,,, = 27V^/3V- y ^[^i + 2Tln(l + e'^^^ + EZ + 2Tln(l + e~f^^)]. (2.54) 
Finally, we obtain the familiar expression of the thermodynamic potential 

+ 2Tln(l + e-^^i) + + 27'ln(l + 6"^^^)]. (2.55) 

The two gap equations A and m in Eqs. (2.32) and (2.36) can be equivalently derived by 
minimizing the thermodynamic potential Eq. (2.55) with respect to m and A, respectively, 
i.e.. 



E. Competition between the chiral and diquark condensates 

In this subsection, through numerical calculations, we will investigate the phase structure 
along the chemical potential direction, analyze the competition mechanism between the 
chiral condensate and diquark condensate, and discuss the influence of the color breaking 
on the quarks with different colors. 

Before the numerical calculations, we should fix the model parameters. The current 
quark mass mo = 5.5MeV, the Fermion momentum cut-off A/ = 0.637GeV, and the cou- 
pling constant in color singlet channel Gs = 5.32GcV~^ arc determined by fitting vac- 
uum properties. The corresponding constituent quark mass in the vacuum is taken to be 
m(/x = 0) = 330MeV. The coupling constant in the color anti-triplet channel Gd can in 
principle be determined by fitting the nucleon properties, here, we choose Gc/Gs as a free 
parameter. 

The double broken phase 

Firstly, we investigate the phase structure along the chemical potential direction with 
respect to different magnitude of Gd/Gs- In the explicit chiral symmetry breaking case, 
we define the point at which the chiral condensate has maximum change as the critical 
chemical potential ji^ for the chiral phase transition, and the point at which the diquark 
condensate starts to appear as the critical chemical potential /ia for the color superconduc- 
tivity phase transition. The two gaps m (white points) and A (black points) determined 
by Eqs. (2.32) and (2.36) are plotted in Fig. 1 as functions of ji with respect to different 
Gd/Gs — 2/3, 1, 1.2, 1.5 in (a), (6), (c) and (d), respectively. 

In Fig. 1 (a), i.e., in the case of Gu/Gs = 2/3, we sec that the chiral phase transition 
and the color superconductivity phase transition occur nearly at the same chemical potential 

~ /^A = 340 MeV. The two phase transitions are of first order. In the explicit chiral 
symmetry breaking case, there is a small chiral condensate in the color superconductivity 
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FIG. 1. The two gaps m (white points) and A (black points) as functions of chemical potential 
H for Gd/Gs = 2/3, 1, 1.2, 1.5, respectively. 



phase /i > fj,^. This phenomena has been called the coexistence of the chiral and diquark 
condensate in Ref. [45]. In this coexistence regime, the chiral condensate is small and can 
be described by the Anderson theorem [102]. 

In Fig. 1 (6), with Gd/Gs — 1, the diquark condensate starts to appear at //a = 298 MeV, 
and the chiral symmetry restores at = 304.8 MeV. Both and fi^ are smaller than 
those in the case of Gd/Gs = 2/3. In the regime from /iA to /x^, we also see the coexistence 
of the chiral and diquark condensates. However, the chiral condensate in this regime is 
not due to the explicit symmetry breaking but due to the dynamical symmetry breaking. 
In order to differ the coexistence of the diquark condensate and the small quark-antiquark 
condensate in the chiral symmetric phase as in Fig. 1 (a), we call the coexistence of the 
diquark condensate and the large quark-antiquark condensate in the chiral symmetry broken 
phase as the "double broken phase" . The diquark gap increases continuously from zero to 
82 MeV in this double broken phase, and jumps up to 152 MeV at the critical point fj,^. 

In Fig. 1 (c) and (d), we see that with increasing of Gd/Gs, the diquark condensate 
starts to appear at smaller fi^, and chiral symmetry restores at smaller fj,^, while the width 
of the regime of double broken phase, /i;^ — /iA, becomes larger. This phenomena has also 
been found in Ref. [101] in the random matrix model. 

However, as shown in Fig. 2 (6), (c) and (d), in the three cases Gd/Gs — 1, 1.2, 1.5, the 
diquark condensate starts to appear even in the vacuum where rib — with rib the baryon 



19 



number density. This, of course, is not physical, because a nonzero diquark condensate 
means the appearance of a nonzero Majorana mass for quarks in vacuum. This would be 
in conflict with the standard definition of the bayron number in vacuum. The appearance 
of the diquark condensate in the vacuum may indicate that, in the NJL model, the diquark 
coupling strength Gd cannot be very large, and it should be smaller than the coupling 
strength in the quark- ant iquark channel. However, it is worth to mention that the 
appearance of the diquark condensate in the vacuum may be a consequence of the lack of 
deconfinement in the NJL model. 
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FIG. 2. The two gaps m (white points) and A (black points) as functions of scaled baron 
density ni,/no for Gd/Gs = 2/3,1,1.2,1.5, respectively. Where uq is the normal nuclear matter 
density. 

We summarize the phase structure along the chemical potential direction: 1) when /i < 
yUA, the chiral symmetry is broken; 2) in the regime from fi^ to fi^, both chiral and color 
symmetries are broken, and 3) when fi > fi^, the chiral symmetry restores partially and 
the phase is dominated by the color superconductivity. The phase structure depends on the 
magnitude of Gd/Gs- If Gd/Gs is very small, the diquark condensate will never appear; 
If Gd/Gs < 1 but not too small, there will be no double broken phase, the chiral phase 
transition and the color superconductivity phase transition occur at the same critical point 
/^x ~ /^A- With increasing of Gd/Gs, the diquark condensate starts to appear at smaller /i^ 
and the chiral symmetry restores at smaller /i^., and the width of the double broken phase, 
f^x~ f^^^ increases. 
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Competition between the chiral and diquark condensates 

In order to understand the competition mechanism and exphcitly show how the diquark 

condensate influences the chiral phase transition, let us analyze Fig. 1 in more detail. In the 
case of Gd/Gs = 0, only chiral phase transition occurs, the thermal system in the mean- 
field approximation is nearly a free Fermi gas made of constituent quarks. In the limit of 
T = 0, there is a very sharp Fermi surface of the constituent quark. When the chemical 
potential is larger than the constituent quark mass in the vacuum, the chiral symmetry 
restores, and the system of constituent quarks becomes a system of current quarks. When 
a diquark gap A forms in the case of Go/Gs 7^ 0, it will smooth the sharp Fermi surface of 
the constituent quark. In other words, the diquark pair lowers the sharp Fermi surface, and 
induces a smaller critical chemical potential of chiral restoration. 

In TABLE I, we fist the chemical potentials //a, at which the diquark gap starts ap- 
pearing, and at which the chiral symmetry restores, for different values of Gd/Gs- 
H% = 345.3 McV is the critical chemical potential in the case of Gd/Gs = 0. A^^ is the value 
of diquark gap at /x^, if there is a jump, it is the lower value. 



Gd/Gs 


/iA(MeV) 


Mx(MeV) 


- /XA(MeV) 


A^(MeV) 




{l4 - l^^)/^x 


1 


298 


304.8 


6 


82 


0.49 


0.57 


1.2 


242 


266 


24 


162 


0.49 


0.64 


1.5 


70 


190 


120 


310 


0.50 


0.89 



TABLE I. The Gd dependence of chemical potentials /xa and //^, jJLp = 345.3MeV. 



We see that for larger Gd/Gs, the diquark condensate appears at a smaller chemical 
potential ^a, and the chiral phase transition occurs at a smaller critical chemical potential 
H^, the gap of the diquark condensate A^ at becomes larger, and the regime of the double 
broken phase becomes wider. 

We assume the relation between A^ and /i^ as 

^^x = I^F- xA^, (2.57) 
and the relation between A^ and //a as 

I^A^I^F- y^x- (2-58) 

In TABEL I, we listed the values of a; = (/Xp — ii^)/Ay^ and y = {^p — /xa)/A;^ for 
different Gd/Gs- It is found that x is almost Gd/Gs independent and equal to 1/2. As for 
y, it is larger than 1/2 and increases with increasing of Gd/Gs- From the Eqs. (2.57) and 
(2.58), we have the relation 

/^x - /^A = (y - 1/2) A^ (2.59) 

for the double broken phase. With increasing of Gd/Gs, y and A^ increase, and then the 
width of the double broken phase, /i^ — /^a, becomes larger. 
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Now we turn to study how the chiral gap influences the color superconductivity phase 
transition. Firstly, we change the constituent quark mass in the vacuum from 330 MeV to 
486MeV. To fit the pion properties, the coupling constant in the quark-antiquark channel 
is correspondingly increased from Gs to 1.2Gs- We plot the diquark gap as a function of /i 
in Fig. 3a for the two vacuum masses and for Gd/Gs — 2/3, 1, 1.2, 1.5. We find that for the 
same Gd/Gs, the diquark gap starts to appear at much larger chemical potential /ia when 
the vacuum mass increases from 330 MeV to 486 MeV. 
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FIG. 3. The influence of chiral gap on the color superconductivity phase transition in the case 
of m{iJ = 0) = 486MeV in (a) and m(/x = 0) = in (6). 

Then we withdraw the quark mass, i.e., taking m — even in the vacuum (this, of 
course, is artificial). We plot the diquark g function of /i in Fig. 3 b for = 0) = 

and 330 MeV and for Gd/Gs = 2/3, 1, 1.2, 1.5. We find that for any Gd/Gs, the diquark 
condensate starts to appear at a much smaller chemical potential for m(/.i = 0) = 
compared with m{iJi. = 0) = 330 MeV. 

From Fig. 3 a and h, we can see that the vacuum quark mass only changes the critical 
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point of color superconductivity //a, the diquark gaps for different vacuum quark masses 
coincide in the overlap regime of the color superconductivity phase, where chiral symmetry 
restores partially. 

Prom the influence of the diquark gap on the chiral phase transition and the influence 
of the chiral gap on the color superconductivity phase transition, it is found that there 
docs exist a strong competition between the two phases. The competition starts at //a and 
ends at Wc call the double broken phase as the competition regime, which becomes 
wider with increasing of Gd/Gs- This competition regime is the result of the diquark gap 
smoothing the sharp Fermi surface of the constituent quark. If the attractive interaction 
in the diquark channel is too small, there will be no diquark pairs, the system will be in 
the chiral breaking phase before fi^, and in the chiral symmetry restoration phase after fi^. 
If the attractive interaction in the diquark channel is strong enough, diquark pairs can be 
formed and smooth the sharp Fermi surface, and induce a smaller critical chemical potential 
of chiral symmetry restoration. 

The influence of color breaking on quark properties 

Finally, wc study how the diquark condensate influences the quark properties. In the 
normal phase, the quarks in different colors are degenerate. However, in the color breaking 
phase, the red and green quarks are involved in the diquark condensate, while the blue 
quarks do not. 
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FIG. 4. The difference of the chiral condensates for quarks in different colors 5 as a function of 
the chemical potential /it with respect to Gd/Gs = 2/3, 1, 1.2, 1.5, respectively. 

The quark mass m appeared in the formulae of this section is the mass for the third 
quark which does not participate in the diquark condensate. We have seen that the diquark 
condensate influences much the quark mass m in the competition regime /ia < fJ^ < fJ^x- 
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the color breaking phase, i.e., when /J, > ha, the quark mass m in different cases of Gd 
decreases slowly with increasing /i, and reaches the same value at about fi = 500MeV. 

The difference of the chiral condensates for quarks in different colors 6 defined in Eq. 
(2.38) is shown in Fig. 4 as a function of the chemical potential fj, with respect to Gd/Gs = 
2/3, 1, 1.2, 1.5. It is found that in any case 5 is zero before /ia, then begins to increase at /ia 
and reaches its maximum at /j,^, and starts to decrease after > /x^, and approaches to zero 
at about /i = 500 MeV, when fi > 500 McV, S becomes negative. With increasing Gd, S{fi^) 
increases from 1 MeV for Gd/Gs = 2/3 to 13 MeV for Gd/Gs = 1.5. Comparing with the 
magnitude of the diquark condensate, 5 is relatively small in the color superconductivity 
phase. 

Summary of this section 

In the mean-field approximation of the extended NJL model, considering only the attrac- 
tive interactions in the 0^ color singlet quark- ant iquark channel and in the color anti-triplet 
diquark channel, the Nambu-Gorkov form of the quark propogator has been evaluated with 
a dynamical quark mass. The Nambu-Gorkov massive propagator makes it possible to 
extend the Green-function approach to the moderate baryon density regime. The familiar 
expression of the thermodynamic potential has been re-evaluated by using the massive quark 
propagator. 

The phase structure along the chemical potential direction has been investigated. The 
system is in the chiral symmetry breaking phase below /xa, in the color superconducting 
phase above /x^, and the two phases compete with each other in the double broken phase 
with width jj,^ — Ha- The existence of the double broken phase depends on the magnitude of 
Gd/Gs- If Gd < Gs, there is no double broken phase, the chiral phase transition and the 
color superconductivity phase transition occur at the same chemical potential. The double 
broken phase starts to appear when Gd/Gs — 1, and the regime of this phase becomes 
wider when the diquark coupling strength increases. 

The double broken phase is a consequence of the competition between the chiral con- 
densate and the diquark condensate. The competition mechanism has been analyzed by 
investigating the influence of the diquark condensate on the sharp Fermi surface of the con- 
stituent quark. The diquark condensate smoothes the sharp Fermi surface, and induces the 
chiral phase transition and the color superconducting phase transition occurring at smaller 
critical chemical potentials. The stronger the diquark coupling strength is, the smaller crit- 
ical chemical potential (or the critical baryon number density) will be for the appearance of 
the diquark condensate. 

However, it is found that when Gd/Gs > 1, the diquark condensate appears even in the 
vacuum. It indicates a nonzero Majorana mass for quarks in the vacuum, which would be in 
conflict with the standard deflnition of the bayron number in the vacuum. This unphysical 
result may suggest that in the NJL model, the diquark coupling strength Gd cannot be 
very large. For the parameters of the NJL model used in this section, our results suggest 
that the diquark coupling strength should be smaller than the coupling strength Gs in the 
quark-antiquark channel. Correspondingly, the upper limit of the diquark gap parameter A 
is about 200 MeV. It is also worth to mention that the appearance of the diquark condensate 
in the vacuum may be a result of the lack of deconflnement in the NJL model. 
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III. THE GAPLESS 2SC PHASE 



In the previous section, we discussed the phase transition from the hadronic phase to 
the color superconducting phase. Now we are going to discuss the two-flavor color super- 
conducting phase under the local charge neutrality condition based on Refs. [87-89,121]. 

It is very likely that the color superconducting phase may exist in the core of compact 
stars, where bulk matter should satisfy the charge neutrality condition. This is because bulk 
matter inside the neutron star is bound by the gravitation force, which is much weaker than 
the electromagnetic and the strong color forces. Any electric charges or color charges will 
forbid the formation of bulk matter. In addition, matter inside neutron star also needs to 
satisfy the /9-equilibrium. 

In the ideal two-flavor color superconducting (2SC) phase, the pairing u and d quarks 
have the same Fermi momenta. Because u quark carries electric charge 2/3, and d quark 
carries electric charge —1/3, it is easy to check that quark matter in the ideal 2SC phase 
is positively charged. To satisfy the electric charge neutrality condition, roughly speaking, 
twice as many d quarks as u quarks are needed. This induces a large difference between 
the Fermi surfaces of the two pairing quarks, i.e., — Ai„ = yUg ~ A*/4, where fi, fie are 
chemical potentials for quarks and electrons, respectively. Naively, one would expect that 
the requirement of the charge neutrality condition will destroy the ud Cooper pairing in the 
2SC phase. 

Indeed, the interest in the charge neutral 2SC phase was stirred by the paper Rcf. [119]. 
It was claimed in this paper that there will be no 2SC phase inside neutron star under the 
requirement of the charge neutrality condition. In fact, the authors meant that for a charge 
neutral three flavor system, the 2SC-|-s phase is not favorable compared to the CFL phase. 
This is a natural result under the assumption of a small strange quark mass, even without 
the requirement of the charge neutrality condition. In the framework of the bag model, in 
which the strange quark mass is very small, the CFL phase is always the ground state for 
cold dense quark matter, and there is no space for the existence of two-flavor quark matter. 

However, there is another scenario about the hadron-quark phase transition in the frame- 
work of the SU(3) NJL model. In the vacuum, quarks obtain their dynamical masses induced 
by the chiral condensate, u, d quarks have constituent mass around 330McV, while the s 
quark has heavier constituent mass, which is around SOOMeV. With the increasing of the 
bayron density, the constituent quark mass starts to decrease when the chemical potential 
becomes larger than its vacuum constituent mass. In this scenario, s quark restores chiral 
symmetry at a larger critical chemical potential than that of u, d quarks. If the deconflne- 
ment phase transition happens sequentially, there will exist some baryon density regime for 
only u, d quark matter and s quark is still too heavy to appear in the system. 

It is worth to mention that the effect of the electric charge neutrality condition on a 
three-flavor quark system is very different from that on a two-flavor quark system. Because 
s quark carries —1/3 electric charges, it is much easier to neutralize the electric charges 
in a three-flavor quark system than that in a two-flavor quark system. However, the color 
charge neutrality condition is nontrivial in a three-flavor quark system, when the strange 
quark mass is not very small. For a detailed consideration of the charge neutral three-flavor 
system, see recent papers Refs. [90-93]. 

In the following, we focus on the charge neutral two flavor quark system. Motivated by 
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the sequential deconfinement scenario, the authors of Ref. [120] investigated charge neutral 
quark matter based on the SU(3) NJL model. To large extent, their results agree with those 
in Ref. [119], i.e., the CFL phase is more favorable than the 2SC+s phase in charge neutral 
three-flavor cold dense quark matter, and they did not find the charge neutral 2SC phase. 

However, it was found in Ref. [121] that a charge neutral two-flavor color superconducting 
(N2SC) phase does exist, which was confirmed in Refs. [122,123]. Comparing with the ideal 
2SC phase, the N2SC phase found in Ref. [121] has a largely reduced diquark gap parameter, 
and the pairing quarks have different number densities. The latter contradicts the paring 
ansatz in Ref. [124]. Therefore, one could suggest that this phase is an unstable Sarma 
state [125]. In Ref. [87], it was shown that the N2SC phase is a stable state under the 
restriction of the charge neutrality condition. As a by-product, which comes out as a very 
important feature, it was found that the quasi-particle spectrum has zero-energy excitation 
in this charge neutral two-flaovr color superconducting phase. Thus this phase was named 
the "gapless 2SC(g2SC)" phase. 

In the following of this section, firstly, using the method introduced in the previous 
section, I will derive the thermodynamic potential for the charge neutral two-fiavor quark 
system based on Ref. [121]. Then I discuss the ground state of the charge neutral two-fiavor 
quark system and introduce the g2SC phase [87]. At last, I show the properties of the g2SC 
phase at zero as well as at nonzero temperatures [88] , and the color screening properties [89] 
in the g2SC phase. 



A. Thermodynamic potential of the neutral two-flavor quark system 

We consider the SU(2) NJL model by assuming that the strange quark does not appear 
in the system, and we only consider scalar, pseudoscalar mesons and scalar diquark. The 
Lagrangian density has the same form as Eq. (2.1). As it was shown in Ref. [45,102], in 
the color superconducting phase, the small quark mass due to the explicit chiral symmetry 
breaking has little effect on the diquark condensate. Therefore, for simplicity, in this section, 
we only consider the chiral limit case mo = 0, and the parameters Gs and A in the chiral 
limit are fixed as 

Gs = 5.0163GeV-^ A = 0.6533GeV. (3.1) 

The corresponding effective mass m — 0.314GeV, and we still choose G^i as a free parameter 
with the ratio r] — Gd/Gs- 

The partition function of the grand canonical ensemble can be evaluated by using stan- 
dard method [114,115], 

Z^N' J[dq][dq]exp{J^ dr J d^x {C + nqjoq)}, (3.2) 

where P — 1/T is the inverse of the temperature T, and is the chemical potential. When 
the electric and color charge neutrality conditions are considered, the chemical potential /i 
is a diagonal 6x6 matrix in flavor and color space, and can be expressed as 

fj, — diag{fJ,ur, l^'ug, l^'ub, IJ'dr, l^'dg, l^'db), (3-3) 
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the chemical potential for each color and flavor quark is specified by its electric and color 
charges 

2 

l^ij,a/3 = il^^ij — l^eQij)^al3 + ^^A'-S^ii (^8)a/3; (3-4) 

where Q and are generators of U{1)q and U{l)s. Here we assumed /X3 = 0, because 
the diquark condenses in the anti-blue direction, and the red and green colored quarks are 
degenerate. For the same flavor, the difference of chemical potentials between the red and 
green quarks and the blue quark is induced by yUg. While for the same color, the difference 
of chemical potentials between u and d is induced by He- 

The explicit expressions for each color and flavor quark chemical potential are given as: 

_ _ 2 1 

_ _ 1 1 

l^dr — IJ'dg — 1^ + g/^e + g/^S) 

2 2 

12 , , 

l^dh = ^l^e- -l^s- (3.5) 

For the convenience of calculations, we define the mean chemical potential Ji for the 
pairing quarks qur, Qdg, and Qug, Qdr 

IJ' = ^ = = - g/^e + gA^S, (3.6) 

and the difference of the chemical potential S/i 

^ fidg - fiur ^ ^dr " ^ug ^ ^3_^^ 

Because the blue quarks do not participate in the diquark condensate, the partition 
function can be written as a product of three parts, 

Z — Zconst^b^rg- (3-8) 

Where the constant part is 

/"^ f (j^ A*A 

Zconst = iV'exp{- dr J ci^x + ^}, (3.9) 

Zb part is for the unpairing blue quarks, and Zrg part is for the quarks participating in 
pairing. In the following, we will derive the contributions of Zf, and Z^g. 

Calculation of Zf, 

Introducing the 8-spinors for quj, and qab, 

= {qub, qdh Qub^ Qdb)^ (3-10) 
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we can express Zi, as 



n,p 

= Det^/2(/3[Go-^],), 



(3.11) 



where the determinantal operation Det is to be carried out over the Dirac, color, flavor and 
the momentum-frequency space, and [Gq has the form of 



( [GtMl 



\ 







\G-,\ll 



\ 






(3.12) 



with 



±1-1 

ia 



(3.13) 



Here we have used ^ = "PpTI^ and = fiialo- 

For the two blue quarks not participating in the diquark condensate, from Eq. (3.11), 
one can have 



1 



ln[Det(/5[G+];,^)Det(/3[Go]„,^)][Det(/3[G+],t)Det(/3[Go],,i)]. 



(3.14) 



Using the results of Eq. (2.42) in the previous section, one has 

[^^mGt]-J)BemG^]-J)]= P'IpI - E^^^U - 
[Det(/?[G+],t)Det(/5[Go-],t)]= P'\pl - E^M - E, 



ub J ' 
'db ] ' 



(3.15) 



with E^^ = E± fi^b and E^^^E± 

lidb where E — ^/p^ + m?. Considering the determinant 
in the flavor, color, spin spaces and the momentum-frequency space, we get the expression 



InZ, = J]{ln(/32bo - {E^b?]) + HP'\pI ' {E^,?]) 

n p 

+\nW, - {E^n + - {EM- 



(3.16) 



Calculation of Z^g 

The calculation of Zrg here is much more complicated than that when the two pairing 
quarks have the same Fermi momonta [51]. Also, we introduce the Nambu-Gokov formalism 
for qur,qug,<ldr and q^g, i.e., 

* = (g«r, qug, qdr, qdg] qur^ qug, Qdr, Qdg)- (3-17) 

The Zrg can have the simple form as 
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where 



with 



n,p 

= Det^/^(/3G;i), 



G 



rg 



i^O ]rg 



A- 



Irg 



rg 






V 










uq 

^^±^-1 



\G^\-^g I 








-,±1-1 

\dg 



and the matrix form for is 



A- = -iA75 



/O 1\ 
0-10 
0-100 
\ 1 0/ 



A+ = 7°(A- 



Prom Eq. (3.18), one obtains 



(3.18) 



(3.19) 



(3.20) 



(3.21) 



InZ,, = -lnDet(/3G'-i). 

Using the identity Eq. (2.45), one can have 

Det(G-^) = DetL>i = Det(-A+A- + A+[Go+]-^[A-]-^[Go ];/) 
= DetiP^ = Det(-A-A+ + A-[Go l^/lA+l-MCol;/)- 

Using the massive energy projectors A± in Eq. (2.21) for each flavor and color 
can work out the diagonal matrix Dy and D2 as 



(3.22) 



(3.23) 
quark, we 



(^i)ii 
(^^1)33 



(/^l)22 = [(PO + 5/i)' - \E-^^\^^ + [(PO + 5/i)' - [i?l]']AH 
(Di)44 = [(PO - 5/x)' - [i?A]']A_ + [(po - <5/x)' - \EV\\^, 



2 11 



2 22 



[{po - 6f,r - [E+]2]A_ + [(po - 5/x)2 - [E^]2]Ah 



where = + 

With the above equations, Eq. (3.22) can be expressed as 

InZ,, = 2 J2 J2i^^i(^"(Po - {El + 5i,)f] + H(3\pI - {E^ - 6^))'] 

n p 

+ HP'ipl - {Ei + Sf,))'] + HP'ipl - {EX - Sf,))']}. 



(3.24) 



(3.25) 
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The thermodynamic potential 

Using the helpful relation Eq. (2.50), one can evaluate the thermodynamic potential of 

the quark system in the mean-field approximation. The total thermodynamic potential for 
u, d quarks in /3-equilibrium with electrons takes the form [121,87,88]: 



127r2 V"^' ^^15 ; 4^5 

+ 3£-E/(gs[^« + 2r,„(l + e-/n], (3.26) 

where the electron mass was taken to be zero, which is sufficient for the purposes of the 
current study. The sum in the second line of Eq. (3.26) runs over all (6 quark and 6 
antiquark) quasi-particles. The explicit dispersion relations and the degeneracy factors of 
the quasi-particles read 

= £;(p) ± [xl] (3.27) 
E%^E{p)±ix^, [xl] (3.28) 

El^^El{p)±^ix. [x2] (3.29) 



B. Gap equations and number densities 

Prom the thermodynamic potential Eq. (3.26), we can derive the gap equations of the 
order parameters m and A for the chiral and color superconducting phase transitions, as 
well as the number densities for quarks with different color and flavor. 

Gap equation for the quark mass 

The gap equation for the quark mass can be derived by using 



9^u,d,e 



dm 

The explicit expression for the above equation is 



= 0. (3.30) 



m = AGs 



d^p m 
{27iyE 



2^[l-/(£^A+)-/(^A-)] 



+2^[i -/(£;+,) -/(£;+_)] 

+2 - fiK,) - He-,) - hE^) - hE-J 



(3.31) 



In the chiral limit, the trivial solution m = corresponds to a chirally symmetric phase of 
quark matter, while a nontrivial solution m ^ corresponds to a phase with spontaneously 
broken chiral symmetry. 

Gap equation for the diqueirk gap 
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Similarly, the gap equation for the diquark condensate can be derived using 

9^u,d,e 



OA 



-0, 



(3.32) 



the explicit expression has the form 



A[l - AGd 



A 



+2-^(1 -/>+,)-/(£;+_))] -0. 



(3.33) 



This equation can also have trivial as well as nontrivial solutions for A, and the nontrivial 
solution A 7^ corresponds to a color superconducting phase. 

Number densities of quarks 

According to one of the criteria of the g2SC phase, the densities of the quark species 
that participate in pairing dynamics are not equal at zero temperature [87,121]. This is in 
contrast to regular pairing in the conventional gapped color superconductors [124] . 

The quark densities are obtained by taking the partial derivatives of the potential flu,d,e 
in Eq. (3.26) with respect to the chemical potentials for each flavor and color. Because of 
the residual SU (2)c symmetry in the ground state of quark matter, the densities of red and 
green quarks for the same flavor, are equal in the g2SC ground state. For example, the 
densities of the up quarks participating in the Cooper pairing read 



n. 



ug 



(27r)3 
E-jl 
E7 



^[l-hEl,)-f{EX_ 
[l-f{E-^)-f{E-_)] 



(3.34) 



Similarly, for the densities of the down quarks, we get 



(27r)3 



E + fi 
Et 



[1 - f{EX.) - f{EX-)] 



^ -^[l-f{E-,)-f{E-_)] 



e; 



-f{E-^) - f{EX^) + hEX-) + f{E-_) 



(3.35) 



As we see, the densities of the up and down quarks participating in the Cooper pairing are 
not equal. In fact, the difference of the densities is given by 



^dg ^ur 



rtdr 

2 



ug 



(27r)3 



f{E-_) - f{E-^) + f{Ei_) - f{Ei^) 



(3.36) 
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which is always nonzero at finite temperature provided the mismatch parameter Sn is 
nonzero. It is even more important for us here that this difi^erence is nonzero at zero tem- 
perature in the g2SC phase of quark matter, 



n. 



urj I J': 



=0 ~ {''^dr '^U3)It=0 



(3.37) 



This is in contrast to the 2SC phase {S/i < A) where this difference is zero in agreement 
with the arguments of Ref. [124]. 

The unpaired bhie quarks, Uf, and d^, are singlet states with respect to SU(2)c symmetry 
of the ground state. The densities of these quarks are 



Hub = 2 

for the up quarks, and 
ridb = 2 



(27r)3 



d^p 



HKb) - fiKb) 



l^ub ( l^ub 

7r2 



(3.38) 



(27r)3 

for the down quarks, respectively. 



HE^b) - hE^b) 



(3.39) 



C. The charge neutral ground state and the g2SC phase 

If a macroscopic chunk of quark matter exists inside compact stars, it must be neutral 
with respect to electric as well as color charges. 

electric and color charge neutrality 

It has been shown that, in QCD case, the two-flavor color superconductivity is automat- 
ically color neutral, because of a non-vanishing expectation value of the 8th gluon field, for 
more details, see Refs. [126,127]. However, in the framework of the NJL model which lacks 
gluons, the color charge neutrality condition is to choose /is such that the system has zero 
net charge ng, i.e.. 



d^p 



-^[^-hE-^.)-f{E, 



Sj (27r)3 
+ ^[l-f{Et,)-f{EX.)] 

+~mb) - hE-b) + ~f{Etb) - ~f{E-j 



A- 



(3.40) 



The detailed numerical calculation in Ref. [121] shows that it is very easy to satisfy the color 
neutrality condition, /xg is around several MeV, which is very small comparing with the quark 
chemical potential ii. In the following discussions, color charge neutrality condition is always 
satisfied. 
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FIG. 5. The effective potential function 
of the electric chemical potential /Xg = 148 MeV 
along the neutrality line (solid line). The results 



of the diquark gap A calculated at a fixed value 
(dashed line), and the effective potential defined 
are plotted for // = 400 MeV with rj = 0.75. 



In the two-flavor quark system, the electric neutrality plays an essential role. Similar to 
the color charge, the electric charge neutrality condition is to choose /le such that the system 
has zero net electric charge uq, i.e.. 



nQ^--ns + p^ + -l^eT^ + 2 ' 



2 ° 37r2 3^^ J (27r)3 

fiE-^) - fiEX^) + f{El_) + f{El 



(3.41) 



To neutralize the electric charge in the homogeneous dense m, d quark matter, roughly speak- 
ing, twice as many d, quarks as u quarks are needed, i.e., ~ 2n„, where Uu^d are the number 
densities for u and d quarks. This induces a mismatch between the Fermi surfaces of pairing 
quarks, i.e., /^d — A*« = A*e = 25//, where //e is the electron chemical potential. 

The proper way to find the charge neutral ground state 

Now, we discuss the role of the electric charge neutrality condition. If a macroscopic 
chunk of quark matter has nonzero net electric charge density ng, the total thermodynamic 
potential for the system should be given by 

^ = ^Coulomb + ^u,d,e, (3.42) 

where ^Couiomb ~ ^qV'^^'^ is the volume of the system) is induced by the repulsive 
Coulomb interaction. The energy density grows with increasing the volume of the sys- 
tem, as a result, it is almost impossible for matter inside stars to remain charged over 
macroscopic distances. So bulk quark matter must satisfy electric neutrality condition with 
^Couiomb\nQ=o = 0, and Qu.d,e\nQ=o is ou the neutrality line. Under the charge neutrality 
condition, the total thermodynamic potential of the system is Q\nQ=o = ^u,d,e\nQ=o- 

Here, we want to emphasize that: The proper way to find the ground state of homogeneous 
neutral u, d quark matter is to minimize the thermodynamic potential along the neutrality 
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FIG. 6. The graphical representation of the solution to the charge neutrality conditions (thick 
dash-dotted line) and the solution to the gap equation for three different values of the diquark 
coupling constant (thick solid and dashed lines). The intersection points represent the solutions 
to both. The thin solid line divides two qualitatively different regimes, A < 6fi and A > dfi. The 
results are plotted for // = 400 MeV and three values of diquark coupling constant Go = i]Gs with 
r] = 0.5, r] = 0.75, and tj = 1.0. 

line Q\nQ=o = ^u,d,e\nQ=o ■ TMs is different from that in the flavor asymmetric quark system, 
where j3- equilibrium is required but fi^ is a free parameter, and the ground state for flavor 
asymmetric quark matter is determined by minimizing the thermodynamic potential ^lu,d,e- 
From Fig. 5, we can see the difference in determining the ground state for a charge 
neutral system and for a flavor asymmetric system when free parameter. In Fig. 5, 

at a given chemical potential /i = 400 MeV and r] = Gd/Gs = 0.75, the thermodynamic 
potential along the charge neutrality line f2|„g=o as a function of the diquark gap A is 
shown by the solid line. The minimum gives the ground state of the neutral system, and the 
corresponding values of the chemical potential and the diquark gap are = 148 MeV and 
A = 68 MeV, respectively. If we switch off the charge neutrality conditions, and consider 
the flavor asymmetric u, d quark matter in (3— equilibrium like in Refs. [128,129], the electric 
chemical potential /ig becomes a free parameter. At a fixed /Xg = 148 MeV and with color 
charge neutrality, the thermodynamic potential is shown as a function of the diquark gap 
by the dashed line in Fig. 5. The minimum gives the ground state of the flavor asymmetric 
system, and the corresponding diquark gap is A = 0, but this state has negative electric 
charge density, and cannot exist in the interior of compact stars. 

Gd dependent charge neutral ground state and the g2SC phase 

Equivalently, the neutral ground state can also be determined by solving the diquark gap 
equation Eq. (3.33) together with the charge neutrahty conditions Eqs. (3.40) and (3.41), 
here we regarded the quark mass is zero in the chiral symmetric phase. 

We visualize this method in Fig. 6, with color neutrality always satisfled, at a given 
chemical potential /i = 400 MeV. The nontrivial solutions to the diquark gap equation as 
functions of the electric chemical potential /Xe are shown by a thick-solid line [t] = Go/Gg = 
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FIG. 7. The quasi-particle dispersion relations at low energies in the 2SC phase (left panel) 
and in the g2SC phase (right panel). 

0.75), a long-dashed line (77 = 1.0), and a short-dashed line (77 = 0.5). It is found that for 
each 77, the solution is divided into two branches by the thin-solid line A = S/i, and the 
solution is very sensitive to r). Also, there is always a trivial solution to the diquark gap 
equation, i.e., A = 0. The solution of the charge neutrality conditions is shown by a thick 
dash-dotted line, which is also divided into two branches by the thin-solid line A — S/J,, but 
the solution of the charge neutrality is independent of rj. 

The cross-point of the solutions to the charge neutrality conditions and the diquark gap 
gives the solution of the system. We find that the neutral ground state is sensitive to the 
coupling constant Gd = V^s in the diquark channel. In the case of a very strong coupling 
(e.g., 77 = 1.0 case), the charge neutrality line crosses the upper branch of the solution to 
the diquark gap, the ground state is a charge neutral regular 2SC phase with A > 611. In 
the case of weak coupling (e.g., 77 = 0.5), the charge neutrality line crosses only the trivial 
solution of the diquark gap, i.e., the ground state is a charge neutral normal quark matter 
with A = 0. 

The regime of intermediate coupling (see, e.g., 7/ = 0.75 case) is most interesting, the 
charge neutrality line crosses the lower branch of the solution of the diquark gap. It should 
be noticed that 77 = 0.75 is from the Fierz transformation, and rj — 2.26/3 ~ 0.75 obtained 
in the SU(2) NJL model in Ref. [106] from fitting the vacuum baryon mass. We will see 
that in the regime of intermediate coupling, the phase with A < 5// is a gapless 2SC (g2SC) 
phase. 

Quasi-peirticle spectrum 

As we already mentioned that the pairing quarks have different number densities in the 
g2SC phase, which is different from the regular 2SC phase. It is the quasi-particle spectrum 
that makes the g2SC phase different from the gapped 2SC phase. 

It is instructive to start with the excitation spectrum in the case of the ideal 2SC phase 
when 5/1 — 0. With the conventional choice of the gap pointing in the anti-blue direction in 
the color space, the blue quarks are not affected by the pairing dynamics, and the other four 
quarsi-particle excitations are linear superpositions of -u^.g and dr^g quarks and holes. The 
quasi-particle is nearly identical with a quark at large momenta and with a hole at small 
momenta. We represent the quasi-particle in the form of Q {quark, hole), then the four 
quasi-particles can be represented explicitly as Q{ur, dg), Q{ug, dr), Q{dr, Ug) and Q{dg, Ur). 
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When Sn = 0, the four quasi-particles are degenerate, and have a common gap A. If there 
is a small mismatch < A) between the Fermi surfaces of the pairing u and d quarks, the 
excitation spectrum will change. For example, we show the excitation spectrum of Q{ur, dg) 
and Q{dg, Uj) in the left panel of Fig. 7. We can see that 5// induces two different dispersion 
relations, the quasi-particle Q{dg,Ur) has a smaller energy gap A — 5/1, and the quasi-particle 
Q{ur, dg) has a larger energy gap A + Sfj,. This is similar to the case when the mismatch is 
induced by the mass difference of the pairing quarks [130]. 

If the mismatch 6fi is larger than the gap parameter A, the lower dispersion relation for 
the quasi-particle Q{dg,Ur) will cross the zero-energy axis, as shown in the right panel of 
Fig. 7. The energy of the quasi-particle Q{dg, Ur) vanishes at two values of momenta p — ii~ 
and p — jJi^ where /j,^ = fl± ^/{SiiY — A^. Thus this phase is called the gapless 2SC (g2SC) 
phase. 

An unstable gapless CFL phase has been found in Ref. [130], and a similar stable gapless 
color superconductivity could also appear in a cold atomic gas [131] or in m, s or s quark 
matter when the number densities are kept fixed [132]. Also, some gapless phases may 
appear due to P-wave interactions in the cold atomic system [133]. 

As we shall see, the interplay of the neutrality condition and the solution to the gap 
equation produces some unusual properties of the g2SC phase that have no analogue in well 
known systems. 



D. The g2SC phase at nonzero temperatures 



In a superconducting system, when one increases the temperature at a given chemical 
potential, thermal motion will eventually break up the quark Cooper pairs. In the weakly 
interacting Bardccn-Copper-Schrieffer (BCS) theory, the transition between the supercon- 
ducting and normal phases is usually of second order. The ratio of the critical temperature 
T^'^^ to the zero temperature value of the gap Aq is a universal value [134] 

TBCS = xscs = — ^ 0-567, (3.43) 



where 7E ~ 0.577 is the Euler constant. In the conventional 2SC phase of quark matter with 
equal densities of the up and down quarks, the ratio of the critical temperature to the zero 
temperature value of the gap is also the same as in the BCS theory [39]. In the spin-0 color 
flavor locked phase as well as in the spin-1 color spin locked phase, on the other hand, this 
ratio is larger than BCS ratio by the factors 2^/^ and 2^/^, respectively. These deviations are 
related directly to the presence of two different types of quasiparticles with nonequal gaps 
[65]. 

This commonly accepted picture of the finite temperature effects in superconducting 
phases changes drastically in the case of dense quark matter when the /3-equilibrium and 
the local neutrality conditions are enforced. Below, we study this in detail. 

Gap equation and charge neutrality condition 

Here we consider the gap equation and charge neutrality conditions in the g2SC phase 
at finite temperature. Because of additional technical complications appearing at finite 
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FIG. 8. The solutions to the gap equation (sohd hues) and the neutrahty condition (dashed 
Unes) calculated for several values of temperature. The results are plotted for fj, = 400 MeV and 
the diquark coupling constant Gd = TjGs with rj = 0.75. In this case, Tc « 50.2 MeV. 

temperature, most of the results that follow will be obtained by numerical computation. It 
is important to keep in mind, however, that our approach remains conceptually the same as 
that at zero temperature. 

Before studying the gap equation (3.33), we recall that the right hand side of this equation 
is a function of the three chemical potentials /x, /le and /ig. The values of /ie and /xg are 
not the free parameters in neutral quark matter. They are determined by satisfying the two 
charge neutrality conditions. In our analysis, on the other hand, it is very convenient to 
enforce only the color charge neutrality by choosing the function A^e? ^) properly. As for 
the electric chemical potential /Ze, it is treated as an independent parameter at intermediate 
stages of calculation. Only at the very end its value is adjusted to make the quark matter 
also electricly neutral. 

As in the case of zero temperature, the resulting color chemical potential /Xg is much 
smaller than the other two chemical potentials in neutral quark matter. This, of course, is 
not surprizing. Small nonzero values of /xg [typically, /xg ~ are required only because 

an induced color charge of the diquark condensate should be compensated. Therefore, 
the smallness of /xg is protected by the smallness of the order parameter. The numerical 
calculations show, in fact, that even the approximation /ig = does not modify considerably 
the exact solutions to the gap equation and to the electric charge neutrahty condition. 
In practice, we either tabulate the function /X8(/i, /le. A) for a given set of parameters, or 
determine it numerically in the vicinity of the ground state solution. 

The solutions to the gap equation (3.33) for several different values of temperature (T = 
0, 20, 40, 50.2 and 60 MeV) are shown graphically in Fig. 8 (solid lines). The values of 
the temperature are marked along the curves. The results are plotted in the (/Xg, A) -plane. 
In computation we kept the quark chemical potential fixed, /i = 400 MeV, and used the 
diquark coupling constant Gd = TjGs with r) = 0.75. This is the same choice that we used 
at zero temperature in Fig. 6. 

Now let us briefly discuss the temperature dependence of the solution. It is not very 
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surprizing that the shape of the graphical solution is smoothed with increasing the temper- 
ature. The same applies to the disappearance of the double-branch structure of the solution 
at a finite value of the temperature, Tg ~ 30 MeV. Indeed, in a model that treats the electric 
chemical potential as a free parameter, this value of the temperature Tg marks the expected 
switch of two regimes. Namely, while the phase transition controlled by the /ig parameter 
(i.e., neutrality is not required) is a first order phase transition at T" < Tg, it becomes a 
second order phase transition at T > T^. 

We also observe that the value of the electric chemical potential at the point where the 
solution to the gap equation intersects with the /Xg-axis has a nonmonotonic dependence on 
the temperature. With increasing temperature, this value increases first and, after reaching 
some maximum value, starts to decrease and goes down to zero eventually. When the 
interplay with the neutrality condition is taken into account later, this simple property of 
the solution would produce rather unusual physical results. 

We note that the temperature dependence of the solutions to the gap equation is very 
similar to the dependence found by Sarma in a sohd state physics analogue of the g2SC 
phase [125]. Of course, there was no analogue of the neutrality condition in the system 
studied in Ref. [125]. Therefore, the mentioned similarity does not extend to the complete 
analysis of the g2SC phase that follows. 

In Fig. 8, we also show the neutrality lines (dashed lines) for the same values of temper- 
ature (T = 0, 20, 40, 50.2 and 60 MeV). The convention is that the lengths of the dashes 
decreases with increasing the value of the temperature. As we see, with increasing temper- 
ature, the neutrality line gets steeper while its intersection point with the /Xg-axis moves 
towards larger values of /ig. The points of intersection of these neutrality lines with the lines 
of solutions to the gap equation, when they exist, are shown as well. 

As in the case of zero temperature, see Fig. 5, we need to show that the points of 
intersections of the gap solutions with the corresponding neutrality lines in Fig. 8 represent 
the ground state of the neutral quark matter. To this end, we calculate the dependence 
of the thermodynamic potential on the vahic of the gap A. Since the charge neutrality 
condition is satisfied only along the neutrality line, we restrict the thermodynamic potential 
only to this line. The numerical results for several values of the temperature are shown in 
Fig. 9. 

The results at ^ = 400 MeV and r] = 0.75 show that, for any T < ^ 50.2 MeV, the 
thermodynamic potential has the global minimum away from the origin, meaning that the 
corresponding ground state develops a nonzero expectation value of A(T). Needless to say 
that this is the same value that one extracts from the geometrical construction in Fig. 8. 
The expectation value A(T) disappears gradually when the temperature approaches Tc from 
below. This is an indication of the second order phase transition. At temperature higher 
than Tc, the ground state of neutral quark matter is the normal phase with A(T) = 0. 

Temperature dependence of the gap 

The temperature dependence of the gap is obtained by numerical solution of the gap 

equation (3.33) together with the two charge neutrality conditions. Of course, this is equiv- 
alent to the geometrical construction used in Fig. 8, where the intersection points of two 
types of lines determine the values of the gaps in the ground state. 
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FIG. 9. The thermodynamic potential of neutral quark matter as a function of the diquark gap 
A calculated for several values of temperature. 

The typical results for the default choice of parameters /x = 400 MeV and rj = 0.75 
are shown in Fig. 10. Both the values of the diquark gap (solid line) and the mismatch 
parameter 6fi = /ie/2 (dashed line) are plotted. One very unusual property of the shown 
temperature dependence of the gap is the nonmonotonic behavior. Only at sufficiently high 
temperatures, the gap is a decreasing function. In the low temperature regime, T < 10 MeV, 
however, it increases with temperature. For comparison, in the same figure, the diquark gap 
in the model with /Xe = and /ig = is also shown (dash-dotted line). This latter has the 
standard BCS shape. 

Another interesting thing regarding the temperature dependences in Fig. 10 appears in 
the intermediate temperature regime, 22.5 ^ T < 37 MeV. By comparing the values of 
A(T) and 6fi in this regime, we see that the g2SC phase is replaced by a "transitional" 2SC 
phase there. Indeed, the energy spectrum of the quasiparticles even at finite temperature is 
determined by the same relations in Eqs. (3.27) and (3.29) that we used at zero temperature. 
When A > 5/1, the modes determined by Eq. (3.29) are gapped. Then, according to our 
standard classification, the ground state is the 2SC phase. 

It is fair to say, of course, that the qualitative difference of the g2SC and 2SC phases is not 
so striking at finite temperature as it is at zero temperature. This difference is particularly 
negligible in the regime of interest where temperatures 22.5 < T < 37 MeV are considerably 
larger than the actual value of the smaller gap, A — 6fi. By increasing the value of the 
coupling constant slightly, however, the transitional 2SC phase can be made much stronger 
and the window of intermediate temperatures can become considerably wider. In either 
case, we find it rather unusual that the g2SC phase of neutral quark matter is replaced by 
a transitional 2SC phase at intermediate temperatures which, at higher temperatures, is 
replaced by the g2SC phase again. 

It appears that the temperature dependence of the diquark gap is very sensitive to the 
choice of the diquark coupling strength rj = Gd/Gs in the model at hand. This is not 
surprising because the solution to the gap equation is very sensitive to this choice. The 
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FIG. 10. The temperature dependence of the diquark gap (sohd hne) and the value of 5fi = /Ue/2 
(dashed hne) in neutral quark matter. For comparison, the diquark gap in the model with = 
and yUg = is also shown (dash-dotted line). The results are plotted for // = 400 MeV and ry = 0.75. 

resulting interplay of the solution for A with the condition of charge neutrality, however, 
is very interesting. This is demonstrated by the plot of the temperature dependence of the 
diquark gap calculated for several values of the diquark coupling in Fig. 11. 

The most amazing are the results for weak coupling. It appears that the gap function 
could have sizable values at finite temperature even if it is exactly zero at zero temperature. 
This possibility comes about only because of the strong influence of the neutrality condition 
on the ground state preference in quark matter. Because of the thermal effects, the positive 
electric charge of the diquark condensate is easier to accommodate at finite temperature. 
This opens a possibility of the Cooper pairing that is forbidden at zero temperature. 

We should mention that somewhat similar results for the temperature dependence of 
the gap were also obtained in Ref. [135] in a study of the asymmetric nuclear matter and 
in Ref. [136] in superfluidity. This behavior can be easily understood: When the cou- 
pling strength is relatively weak, at zero temperature, it is very difficult to form a pair for 
two fermions at largely separating sharp Fermi surfaces. However, when the termperature 
increases, the thermal excitations smooth the Fermi surfaces, thus offer the possibility for 
forming Cooper pairs. Of course, when the temperature is very high, the thermal excitations 
will eventually destroy the Cooper pairs. 

Nonuniversal ratio Tc/ Aq 

In the preceding subsection, we saw that the temperature dependence of the gap in the 
g2SC phase is very different from the standard benchmark result in the BCS theory. One of 
the main properties of the BCS temperature dependence is a universal value of the ratio of 
the critical temperature to the value of the gap at zero temperature Aq, see Eq. (3.43). It 
is instructive, therefore, to calculate the same quantity in the g2SC phase of neutral quark 
matter. 



Let us start from a simple exercise, and consider the results plotted in Fig. 10 more 
carefully. First of all, we find that T^^^ ^ 62.06 MeV and Ag^^ = 109.4 MeV when there is 
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FIG. 11. The temperature dependence of the diquark gap in neutral quark matter calculated 
for several values of the diquark coupling strength rj = Gd/Gs- 

no mismatch parameters in the model (i.e., He = and /xg = represented by dash-dotted 
line). It is not very surprizing that the ratio of interest r2sc ~ 0.567 is in agreement with the 
BCS result. Now, if we check the results for the gap function in the g2SC phase (solid line), 
we find that Tf^^ ^ 50.2 MeV and Af^^ ^ 68.2 MeV. Thus, the ratio is r^asc ^ 0.7357 
which is considerably larger than the 2SC result. 

It appears that the real situation is even more interesting. The result for the ratio 
rg2sc ~ 0.7357 is not universal. In fact, its value depends very much on the diquark 
coupling constant and, moreover, it can even be arbitrarily large. This last statement may 
not be so unexpected if we recall the temperature dependences of the gap shown in Fig. 11. 

The numerical results for the ratio of the critical temperature to the zero temperature 
gap in the g2SC function of the diquark coupling strength t] = Gd/Gs are plotted 

in Fig. 12. The dependence is shown for the most interesting range of values of ?7 = Gd/Gs, 
0.68 < < 0.81, which allows the g2SC stable ground state at zero temperature. When 
the coupling gets weaker in this range, the zero temperature gap vanishes gradually. As we 
saw from Fig. 11, however, this does not mean that the critical temperature vanishes too. 
Therefore, the ratio of a finite value of Tc to the vanishing value of the gap can become 
arbitrarily large. In fact, it remains strictly infinite for a range of couplings. 

E. Chromomagnetic instability in the g2SC phase 

The g2SC phase has four gapless modes and two gapped modes, one may think that the 
low energy (large distance scale) properties in the g2SC phase should interpolate between 
those in the normal phase and those in the 2SC phase. However, its color screening properties 
do not fit this picture. 

As we know, one of the most important properties of the ordinary superconductor is its 
Meissner effect, i.e., a superconductor expels the magnetic field, which was discovered by 
Meissner and Ochsenfeld in 1933 [35]. From theoretical point of view, the Messiner effect 



41 



0.7 0.74 0.78 

11 

FIG. 12. The ratio of the critical temperature to the zero temperature gap in neutral quark 
matter as a function of the coupling strength r] = Gd/Gs- 

can be explained using the linear response theory. The induced current j™'' is related to the 
magnetic field Aj by j*"'^ = IlijA\ where the response function Uij is the photon polarization 
tensor. The response function has two components, diamagnetic and paramagnetic part 
[137]. In the static and long- wavelength limit, for the normal metal, the paramagnetic 
component cancels exactly the diamagnetic component. While in the superconducting phase, 
the paramagnetic component is quenched by the energy gap and producing a net diamagnetic 
response. Thus the ordinary superconductor is a prefect diamagnet. 

In color superconducting phases, the gluon self-energy (the response function to an ex- 
ternal color field), has been investigated in the ideal 2SC phase [54] and in the CFL phase 
[59]. The results show that the gauge bosons connected with the broken generators obtain 
masses in these phases, which indicate the Meissner screening effect in these phases. 

It is very interesting to know the chromomagnetic property in the g2SC phase. We 
studied the g2SC phase in the framework of the SU(2) NJL model, and the NJL model 
lacks gluons. As reflection of this, it possesses the global instead of gauged color symmetry. 
In addition, there appear five Nambu-Goldstone (NG) bosons in the ground state of the 
model when the color symmetry is broken. In QCD, there is no room for such NG bosons. 
However, the NJL model can be thought of as the low energy theory of QCD in which 
the gluons, as independent degrees of freedom, are integrated out. The gluons could be 
reintroduced back by gauging the color symmetry in the Lagrangian density of the NJL 
model, providing a semirigorous framework for studying the effect of the Cooper pairing on 
the physical properties of gluons. 

The existence of the g2SC phase can be regarded as a physical and model independent 
result under the restriction of local charge neutrality condition, the order parameter for this 
phase is A < 6fi. In Ref. [89], we calculated the gluon self-energy in the g2SC phase. It 
is found that, in this phase, the symmetry broken gauge bosons have imaginary Meissner 
screening masses, which is induced by the dominant paramagnetic contribution to the gluon 
self-energy. In condensed matter, this phenomenon is called the paramagnetic Meissner 
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effect (PME) [138], and fias been observed in some fiigfi temperature superconductors and 
small superconductors. 

Unavoidably, the imaginary Meissner screening mass indicates a chromomagnetic insta- 
bility of the g2SC phase. There are several possibilities to resolve this instabihty. One is 
through a gluon condensate to stabilize the system, which may not change the structure of 
the g2SC phase. It is also possible that the instability drives a new stable ground state, 
which may have a rotational symmetry breaking like in Refs. [139,140], or even have an 
inhomogeneous phase structure, like a crystal [68-70] or a vortex [141] structure. In the 
future, one has to resolve this problem. 

IV. THE MIXED PHASE AND NONSTRANGE HYBRID STAR 

We have discussed the homogeneous two-flavor quark matter when charge neutrality 
conditions are satisfied locally, and found that the local charge neutrality conditions impose 
very strong constraints on determining the ground state of the system. In this section, we 
are going to discuss the mixed phase and the nonstrange hybrid star based on Ref. [142]. 

A. The mixed phase 

When charge neutrality conditions are satisfied globally, one can construct a mixed phase 
[142-144]. Inside mixed phases, the charge neutrality is satisfied "on average" rather than 
locally. This means that different components of mixed phases may have nonzero densities 
of conserved charges, but the total charge of all components still vanishes. In this case, one 
says that the local charge neutrality condition is replaced by a global one. There are three 
possible components: (i) normal phase, (ii) 2SC phase, and (iii) g2SC phase. 

The pressure of the main three phases of two-flavor quark matter as a function of the 
baryon and electric chemical potentials is shown in Fig. 13 at r] = 0.75. In this figure, we 
also show the pressure of the neutral normal quark and gapless 2SC phases (two dark solid 
lines). The surface of the g2SC phase extends only over a finite range of the values of /Xg. 
It merges with the pressure surfaces of the normal quark phase (on the left) and with the 
ordinary 2SC phase (on the right). 

It is interesting to notice that the three pressure surfaces in Fig. 13 form a characteristic 
swallowtail structure. As one could see, the appearance of this structure is directly related to 
the fact that the phase transition between color superconducting and normal quark matter, 
which is driven by changing parameter /le, is of first order. In fact, one should expect 
the appearance of a similar swallowtail structure also in a self-consistent description of the 
hadron-quark phase transition. Such a description, however, is possible. 

From Fig. 13, one could see that the surfaces of normal and 2SC quark phases intersect 
along a common line. This means that the two phases have the same pressure along this 
line, and therefore could potentially co-exist. Moreover, as is easy to check, normal quark 
matter is negatively charged, while 2SC quark matter is positively charged on this line. 
This observation suggests that the appearance of the corresponding mixed phase is almost 
inevitable. 



43 





FIG. 13. At r] = 0.75, pressure as a function of = Hb/'^ £ind fi^ for the normal and color 
superconducting quark phases. The dark solid lines represent two locally neutral phases: (i) the 
neutral normal quark phase on the left, and (ii) the neutral gapless 2SC phase on the right. The 
appearance of the swallowtail structure is related to the first order type of the phase transition in 
quark matter. 

Let us start by giving a brief introduction into the general method of constructing mixed 
phases by imposing the Gibbs conditions of equihbrium [145,146]. From the physical point 
of view, the Gibbs conditions enforce the mechanical as well as chemical equilibrium between 
different components of a mixed phase. This is achieved by requiring that the pressure of 
different components inside the mixed phase are equal, and that the chemical potentials (/x 
and fie) are the same across the whole mixed phase. For example, in relation to the mixed 
phase of normal and 2SC quark matter, these conditions read 

P(^«)(/.,/.e)=P(^^^H/i,/ie), (4.1) 
= f,i^Q) = f,i^SC)^ (4 2) 

/^e = /.r^ = /^r^^ (4.3) 

It is easy to visualize these conditions by plotting the pressure as a function of chemical 
potentials (/i and fie) for both components of the mixed phase. This is shown in Fig. 14. As 
should be clear, the above Gibbs conditions are automatically satisfied along the intersection 
line of two pressure surfaces (dark solid line in Fig. 14). 

Different components of the mixed phase occupy different volumes of space. To describe 
this quantitatively, we introduce the volume fraction of normal quark matter as follows: 
X2SC = Vnq/V (notation Xb Cleans volume fraction of phase A in a mixture with phase 
B). Then, the volume fraction of the 2SC phase is given by xj^Q = (1 ~ X2sc)- From the 
definition, it is clear that < X2sc — ^■ 

The average electric charge density of the mixed phase is determined by the charge 
densities of its components taken in the proportion of the corresponding volume fractions. 
Thus, 
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FIG. 14. Ai r] = 0.75, pressure as a function of = and fig for the normal and color 

superconducting quark phases (the same as in Fig. 13, but from a different viewpoint). The dark 
solid line represents the mixed phase of negatively charged normal quark matter and positively 
charged 2SC matter. 



ni^^^ = xS^r^/^, /^e) + (1 - /^e). (4.4) 

If the charge densities of the two components have opposite signs, one can impose the global 
charge neutrality condition, ni^^^^ = 0. Otherwise, a neutral mixed phase could not exist. 
In the case of quark matter, the charge density of the normal quark phase is negative, while 
the charge density of the 2SC phase is positive along the line of the Gibbs construction (dark 
solid line in Fig. 14). Therefore, a neutral mixed phase may exist. The volume fractions of 
its components are 

(2SC) 

NQ _ s 

X2SC - (2SC) {NQ) ' y^-^) 
Tie TT'e 

Xnq = 1 - X2SC = -JWq) (2SC0 ■ '^^■^^ 

Tie — Tie 

After the volume fractions have been determined from the condition of the global charge 
neutrality, we could also calculate the energy density of the corresponding mixed phase, 

.(^^) = xS^(^«)(/^,/^e) + (1 - X^i)e'''''\f^,^^e). (4.7) 

This is essentially all that we need in order to construct the equation of state of the mixed 
phase. 

So far, we were neglecting the effects of the Coulomb forces and the surface tension 
between different components of the mixed phase. In a real system, however, these are 
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important. In particular, the balance between the Coulomb forces and the surface tension 
determines the size and geometry of different components inside the mixed phase. 

In our case, nearly equal volume fractions of the two quark phases are likely to form 
alternating layers (slabs) of matter. The energy cost per unit volume to produce such layers 
scales as a'^/^{nf^^'' — ni^^^y^^ where a is the surface tension [147]. Therefore, the quark 
mixed phase is a favorable phase of matter only if the surface tension is not too large. Our 
simple estimates show that (imax < 20 MeV/fm^. However, even for slightly larger values, 
20 < a < 50 MeV/fm^, the mixed phase is still possible, but its first appearance would occur 
at larger densities, 3po < Ps < 5po- The value of the maximum surface tension obtained 
here is comparable to the estimate in the case of the hadronic-CFL mixed phase obtained 
in Ref. [148]. The thickness of the layers scales as a^/^{nf^^'^ - ni^^^^/^ [147], and its 
typical value is of order 10 fm in the quark mixed phase. This is similar to the estimates in 
various hadron-quark and hadron-hadron mixed phases [147,148]. While the actual value of 
the surface tension in quark matter is not known, in this study we assume that it is not very 
large. Otherwise, the homogeneous gapless 2SC phase should be the most favorable phase 
of nonstrange quark matter [87]. 

Under the assumptions that the effect of Coulomb forces and the surface tension is small, 
the mixed phase of normal and 2SC quark matter is the most favorable neutral phase of 
matter in the model at hand with rj = 0.75. This should be clear from observing the pressure 
surfaces in Figs. 13 and 14. For a given value of the baryon chemical potential /x = /xb/3, 
the mixed phase is more favorable than the gapless 2SC phase, while the gapless 2SC phase 
is more favorable than the neutral normal quark phase. 

B. Nonstrange hybrid star 

In order to construct a neutron star, we need to know the equation of state from very 
low baryon density to very high baryon density. Unfortunately, currently there is no single 
model that can describe the whole bayron density regime very successfully. In the following, 
we use a QCD-motivated hadronic chiral SU(3)i;,xSU(3)ij model in normal baryon density 
regime, and use the NJL model in the quark phase. 

Hadronic matter 

At densities around normal nuclear matter density po ~ 0.15 fm~^, the description of 
baryonic matter in terms of quarks could hardly be adequate. At such low densities, quarks 
are confined inside hadrons. Thus, it is more natural to use an effective hadronic model. 

We use a QCD-motivated hadronic chiral SU(3)/,xSU(3)r model as an effective theory of 
strong interactions to describe the low density regime of the baryonic matter [149-151]. This 
model was found to describe reasonably well the hadronic masses of various SU(3) multiplets, 
finite nuclei, hypernuclei, excited nuclear matter and neutron star properties [149-151]. 
The basic features of the chiral model are: (i) Lagrangian of the model is constructed in 
accordance with the nonlinear realization of the chiral SU(3)iXSU(3)i^ symmetry; (ii) heavy 
baryons and mesons get their masses as a result of spontaneous symmetry breaking; (iii) 
the masses of pseudoscalar mesons (pseudo-Nambu-Goldstone bosons) result from an explicit 
symmetry breaking; (iv) a QCD-motivated field that describes the gluon condensate (dilaton 
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FIG. 15. Pressure as a function of ^ = /is/3 and /Ug for the hadronic phase (at the bottom), 
for the two-flavor color superconducting phase (on the right in front) and the normal phase of 
quark matter (on the left, and back on the right). The dark thick line follows the neutrality line in 
hadronic matter, and two mixed phases: (i) the mixed phase of hadronic and normal quark matter; 
and (ii) the mixed phase of normal and color superconducting quark matter. 

field) is introduced. 

The mixed phase of hadronic matter and quark matter 

It is expected that the phase transition between the hadronic phase and the normal 
quark phase is a first order phase transition at zero temperature and finite baryon chemical 
potential. Then, the hadronic and quark phases could co-exist in a mixed phase [145,146]. 
This mixed phase should satisfy the Gibbs conditions of equilibrium which are similar to 
those discussed in the previous section, see Eqs. (4.1)-(4.3). The total energy density in the 
hadron-quark mixed phase is given by 

,(MP) ^ ;,^«e(^«)(^,/..) + (1 - (4.8) 

where Xh^ denotes the volume fraction of normal quark matter inside a mixture with 
hadronic matter. 

To visualize the Gibbs construction of the hadron-quark mixed phase, we plot the 
hadronic surface of the pressure along with the quark surfaces discussed in the previous 
section. Thus, Fig. 14 is replaced by Fig. 15. The new figure shows the surfaces of the 
pressure of the pure hadronic and quark phases as a function of their chemical potentials 
/i and lie- The intersection lines of different surfaces indicate all potentially viable mixed 
phase constructions. Although the Gibbs conditions are satisfied along all these lines, not 
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all of them could produce globally neutral phases (e.g., there are no neutral constructions 
along the light shaded solid lines in Fig. 15). 

The dark solid line gives the complete, most favorable construction of globally neutral 
matter in general /^-equilibrium. This line consists of three pieces. The lowest part of 
the curve (up to the point denoted by □ and P < 10 MeV/fm^) corresponds to the pure 
confined hadronic phase. Within this regime matter is mostly composed of neutrons with 
little fractions of protons and electrons to realize the charge neutrality and /^-equilibrium. 
The hyperonic particles (A, E and S) are not present in this lowest density regime. Such 
particles would appear in the hadronic phase at considerably higher values of pressure and 
density. 

The mixed phase of hadronic and normal quark matter starts at the baryonic density 
Pb ~ 1.49po which corresponds to the D-point in Fig. 15. At this point the first bubbles of 
deconfined quark matter appear in the system. At the beginning of this hadron-quark mixed 
phase, the deconfined bubbles are small but highly negatively charged, whereas hadronic 
matter, in which the bubbles are embedded, is shghtly positively charged. The global charge 
neutrality condition reads 

- X^V^'^H/^,/^e) + (1 - X^^)nf H/^,/^e) = 0. (4.9) 

where rii^^ and rii^^^ are the charge densities of hadronic and normal quark matter, respec- 
tively. This condition should be satisfied at each point along the middle part of the dark 
solid line (i.e., between the points denoted by □ and A). With increasing density (from 
about 1.49po up to 2.56po), the volume fraction of hadronic matter decreases (down to about 
0.59), while the fraction of normal quark matter increases (up to about 0.41). 

Equation of state 

The equations of state for quark and hybrid matter are shown in Fig. 16. The first 
equation of state corresponds to globally neutral quark matter which is a mixture of the 
normal quark and 2SC phases. As for the equation of state of hybrid matter, it is constructed 
out of the equation of state of neutral hadronic matter and two Gibbs constructions in 
accordance with Fig. 15. As before, the points that indicate the beginning of two different 
mixed phases arc denoted by □ and A in Fig. 16. 

Mass-radius relation for hybrid star 

The mass-radius relations for hybrid and pure quark stars are shown in Fig. 17. As one 
would expect, the pure quark stars composed of the mixed phase have much smaller radii 
and the value of their maximum mass is slightly smaller, see Fig. 17. Our results for pure 
quark stars are comparable to those in Refs. [122,152]. The difference between hybrid and 
pure quark stars is mostly due to the low density part of the equation of state. This is also 
evident from the qualitative difference in the dependence of the radius as a function of mass 
for the hybrid and quark stars with low masses. The corresponding hybrid stars are large 
because of a sizable low density hadronic layer, while the quark stars are small because that 
have no such layers. 

Here we use a two-fiavor version of the NJL model to describe the deconfined phase. 
However, the three-flavor extension of the NJL model in Ref. [120] suggests that strange 



48 




FIG. 16. The equation of state for globally neutral hybrid matter (solid line) and globally 
neutral quark matter (dashed line). The points of the beginning of the two mixed phases are 
denoted by □ and A. 
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FIG. 17. The mass-radius relations for hybrid stars (solid line) and quark stars (dashed line). 
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quarks might be present in matter above a critical density of about 5po- For the consideration 
of strange quark star in the framework of the NJL model, please see Ref. [153]. 

V. CONCLUSION AND OUTLOOK 

This article focuses on the two-flavor color superconducting phase at moderate baryon 
density. 

In order to simultaneously investigate the chiral phase transition and the color super- 
conducting phase transition, the Nambu-Gorkov formalism has been extended to treat the 
quark-antiquark and diquark condensates on an equal footing, and the competition between 
the chiral and diquark condensates has been investigated. 

Dense u, d quark matter under local and global charge neutrality conditions in /3- 
equilibrium has been discussed in detail. Under the local charge neutrality condition, the 
homogeneous neutral ground state is sensitive to the coupling constant in the diquark chan- 
nel, it will be in the regular 2SC phase when the diquark coupling is strong, in the normal 
phase when the diquark coupling is weak, and in the g2SC phase in the case of intermediate 
diquark coupling. The low energy quasi-particle spectrum of the g2SC phase contains four 
gapless modes and only two gapped modes. In the intermediate diquark coupling regime, 
under global charge neutrality condition, assuming that the effect of Coulomb forces and the 
surface tension is small, one can construct a mixed phase composed of a positive charged 
2SC phase and a negative charged normal quark phase. It will be valuable to compare the 
homogeneous g2SC phase and the mixed phase, and give a clear answer to which phase is 
preferable. 

Even though the regime of the diquark coupling strength for the existence of the g2SC 
phase is rather narrow, i.e., 0.7 ^ fj < 0.8, it is noticed that the value rj = 0.75 from 
Fierz transformation and also the value rj ~ 2.26/3 from fitting the baryon vacuum mass, 
are inside this regime. The g2SC has rather unusual properties at zero as well as at finite 
temperature, and its chromomagnetic properties are extremely surprising. The paramagnetic 
Meissner effect indicates a chromomagnetic instability in the g2SC phase. In the near 
future, we need to resolve the problem of the chromomagnetic instability in the g2SC phase, 
and to investigate whether the chromomagnetic instability will drive a new inhomogeneous 
structure. 

In this article, I only discussed the two-flavor quark system. In the moderate baryon 

density regime // ~ 500MeV, the strange quark may appear in the system. The cold dense 
three-flavor quark system has a more complicated and more interesting phase structure. 
In different conditions, the three-flavor quark system can be in one of the phases of CFL, 
CFL+K, CFL+r/, gCFL, and dSC [78,79,86,90-94]. It is, of course, very interesting and 
also challenging to investigate all these phases in one self-consistent framework. For the 
charge neutral gCFL phase, it is very attractive to know whether there is chromomagnetic 
instability. 

In our calculations, we used the mean-fleld approximation. It is very interesting to see 
how the diquark fluctuation will affect the results. It should be noticed that till now, we are 
using the BCS scenario even in the intermediate baryon density regime. From the experience 
in condensed matter, the BCS scenario only works when the attractive interaction is very 
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weak, where the Cooper pairs are loosely bound states. While with the increasing of the cou- 
pling strength, the Cooper pairs become tightly bounded, the fluctuation from the Cooper 
pairs cannot be neglected. In the strong coupling case, the Cooper pairs become real bound 
sates in the coordinates space, and a Fermion system transfers to a boson system. The 
boson system may experience the Bose- Einstein condensation (BEC). This is the basic idea 
of the crossover from BCS to BEC [154,155]. The fluctuation of the gauge flclds as well as 
of the diquark were investigated in Refs. [156,157], and it will be very interesting to see how 
the Cooper pair structure [158] changes when fluctuations are considered. In addition to the 
diquark fluctuation, the quark- antiquark fluctuation may also need to be considered in the 
intermediate baryon density regime. It was found in Ref. [159] that the quark-antiquark fluc- 
tuation has signiflcant influence on the chiral phase transition. When the quark-antiquark 
fluctuation is considered self-consistently [160], the critical baryon density for the chiral 
phase transition becomes much larger than that in the mean-fleld approximation. 

The most interesting thing in the held of color superconductivity, of course, is to look 
for this phase in the nature or to create this phase in the laboratory. There are already 
some work to flnd the possibility of the existence of the color superconducting phase inside 
the neutron star, like the mass-radius relation [122,142,152,153] and the cooling history 
[161]. The energy release from a hadronic star to a color superconducting star may relate 
to the supernovae or the gamma-ray burst [162]. About the possibility of creation color 
superconductivity in heavy-ion collisions, the work is in progress [163]. Nevertheless, we are 
just at the very beginning of this fleld, and there is still a long way to go. 
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APPENDIX A: FIERZ TRANSFORMATION 



Color-current interaction has the form of 

8 



(Al) 



0=1 



where = (A")c/2 are the SU(3) color generators with tr{tH^) = bahj'^-, a, /5, 7,(5 are color 
indices acting on color space, i,j,k,l are flavor indices acting on flavor space, and r,s,t,u 
are spinor indices acting on Dirac space. 

Fierz transformation in the Dirac space 

Defining 5 = 1 (g) 1, y = 7^ (g) 7^^, T = cr^^ ® af"", A = 7^75 (g) 7'^75, P = 75 ® 75, the 
exchange of the spinor indices s and or r and t, can be represented as 



^ rs,tu / J ^%J^ ru,tsi 



(A2) 



where V' e {S, V, T, A, P), and 



CijCjk — Sih- 



(A3) 



The explicit expression for the above equation reads 



/-1/4 -1/4 -1/8 1/4 -1/4 \ 

-1 1/2 1/2 1 

-3 1/2 -3 

1 1/2 1/2 -1 



/ s\ 


1 


V 




T 




A 









( s\ 

V 

T 
A 



\-l/4 1/4 -1/8-1/4-1/4/ 
In the following, we only need the transformation of V , i.e., 

® l^)rs,tu = -(1 «) 1 + i75 ® ilb)ru,ts 



(A4) 



(A5) 



where = {1, ^75, -^Ifi, ■^Ifi'lb}, from which we can see that the contribution of the tensor 
part is zero, and thus the chiral symmetry is preserved. 
For diquarks, we have 

= Cy,aCu{ln)rs{l'^)ab 

^-{KaC®CK%t,^s, (A6) 



where C — i7^7° is the charge conjugation matrix. 
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1. Fierz transformation in the quark-antiquark sector 

We will use the Fierz transformation in the color, flavor and Dirac spaces respectively. 
Fierz identity in the color space 
The standard identity 

c c 1 1 1 c c 

VlS " 2^"^ ~ Jp)SasS'yl3 - Jj-t'ast^^l3, (A7) 

gives both the color singlet and color octet terms. 
Fierz identity in the flavor space 

Rewrite 

e 

where has been ordered in the way flavor singlet + octet, with 
Color singlet psirt 

Using Eqs. (A5), (A7) and (A8), we can get the color singlet part of the Lagrange Eq. 
(Al), 

a=0 

- ^G,^[(^7.f V')^ + (H.ls^m- (AlO) 

a=0 



2. Fierz transformation in the diquark sector 



Fierz identity in the color space 

taptjs = XI Ola ~ ^t-ytip^ (All) 

S A 

where 5" e {0, 1, 3, 4, 6, 8} denotes the symmetric and A e {2, 5, 7} the antisymmetric Gell- 

Mann matrices. 

The above identity is decomposed in the way of 6c + 3c, which shows that the interaction 
in the 6c channel is oppositive to that in the 3c channel. There is another decomposition in 
the way of Ic + 3c, i.e.. 
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= ^(1 - Yj^asS.p + ^/a,4p, (A12) 



which shows that the interaction in the color antitriplet channel has the same symbol as 
that in the singlet channel. The interaction in the color singlet channel Ic is attractive, so 
is the color antitriplet channel Sc- The former is responsible for the bound state of meson, 
the latter is responsible for the bound state of baryon, as well as the quark-quark Cooper 
pairing at high baryon density. The ratio of the magnitude of the coupling constant in 3c 
channel and Ic channel is 1/2 from Eq. (A12). 

Fierz identity in the flavor space 

Sij^ki = Y,iG')ik{G')ij, (A13) 



where 



X2 X5 X7 1 \0 \1 \3 \4 \6 \8 

lr<c\ _ jr^e r^ei _ r Af _ Af Af Af Af ■^F -j 

^ i - - 1^, ^, ^1 - ^1- 



(A14) 



where -S" e {0, 1, 3, 4, 6, 8} denotes the symmetric and A e {2, 5, 7} the antisymmetric Gell- 
Mann matrices. 

Color anti-triplet part 

First, using Eqs. (A6), we have 

= [^|;KaC^|J'^][VCK^^|^]. (A15) 
Then, using Eq. (A12) the color anti-triplet part can be expressed as 

['^|Jar{7t^)rs{t'')ap^Ps[i>^tiY)tu{t"')J^i>Su\ 

= —m^c'^ni^^CK^'-^i^]. (A16) 

At last, the flavor structure is determined by the Pauli principle. A color triplet diquark 
state is antisymmetric under the exchange of colors, due to the Pauli principle, the matrices 
CK"-\f must be totally symmetric. Thus the axial- vector and scalar diquark currents trans- 
form according to antisymmetric flavor 3 representation, whereas the scalar-, pseudoscalar- 
and vector-diquarks are necessarily flavor sextetes. 



APPENDIX B: CHARGE CONJUGATION 

Charge conjugation is conventionally deflned to take a fermion with a given spin orien- 
tation into an antifermion with the same spin orientation. 
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The charge-conjugated filed ^'^ is defined by 

V^(x) = Ci:^^{x),i:{x) = i^^^{x)C, (Bl) 
where C — i^'^^^ is the charge-conjugation matrix in the Dirac representation, 

C = -C-^ = -C^ = -C^ (B2) 

and 

C7^C-' = -{-i^f. (B3) 

The charge conjugation of bihnears is a httle tricky, and it helps to write in spinor indices. 
For the scalar, 

^^^^ - € 

= -0^- (B4) 

In the last step, the fermion anticommunication and CC — —1 have been used. 
For the vector, 

= -tpjf'ip. (B5) 



For the kinetic term. 



^V^M^"" = (B6) 



usmg 



d.ii^r^) = id,i,)r^ + (B7) 
and the vector current conservation, we have 

ijj'^Yd^i^^ = ^7^5/.^- (B8) 
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